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SECTION  1 
INTRODUCTION 


! 

,1 


There  are  a  wide  variety  of  electro-optical  sensors  employed  in 
the  modem  tactical  battlefield  environment.  The  environment  in  which 
the  sensors  must  operate  is  often  extremely  adverse.  The  two  categories 
of  degradation  sources  which  can  significantly  affect  sensor  performance 
are  the  natural  and  the  man-made.  The  natural  sources  include  the  ad¬ 
verse  weather  effects  such  as  rain  and  fog,  the  effects  of  naturally 
occurring  aerosols  such  as  wind  driven  dust,  and  the  natural  thermal 
environment.  The  man-made  sources  include  all  anthropogenic  gases, 
particles,  and  thermal  effects  produced  by  human  activity.  These  in¬ 
clude  countermeasures  smoke,  vehicular  exhaust  and  dust,  muzzle  blast 
emissions,  explosion  produced  products  such  as  gases,  combustion  par¬ 
ticles,  and  thrown-up  dust  and  vegetation,  and  smoke  from  burning 
vehicles,  structures,  and  vegetation.  In  addition  to  the  gases  and 
particles  introduced,  explosions  and  fires  are  a  source  of  heat  and 
atmospheric  convection.  The  resulting  thermal  environment  can  severely 
degrade  sensor  performance. 

The  propagation  effects  due  to  the  natural  environment  have  been 
studied  and  modeled  for  many  years;  in  the  last  few  years  most  of  the 
man-made  effects  have  received  extensive  attention.  Experimental  test 
programs  have  been  run  for  countermeasures  smoke  and  explosion  produced 
dust  clouds.  Models  have  been  or  are  in  the  process  of  being  developed 
for  most  effects.  In  this  report  we  begin  the  study  of  the  thermal 
effects  of  battlefield  fire  plumes. 

Fires  in  a  battlefield  can  be  localized,  such  as  burning  vehicles 
and  structures,  or  area  fires,  such  as  a  burning  field.  Both  the  active 
flame  and  the  rising  heated  air  plume  above  the  flames  are  regions  which 
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can  produce  signal  propagation  degradation.  Here  we  are  concerned  with 
the  rising  heated  air  plume.  Turbulence  in  the  hot  rising  plume  can 
cause  signal  amplitude  and  phase  fluctuations,  defocusing,  beam  wander¬ 
ing  (jitter)  and  spreading,  time  delay,  and  pulse  distortion.  In  add: 
tion  the  gases  and  particulate  combustion  products  in  the  plume  can 
cause  absorption  and  scattering.  In  this  report  we  derive  estimates 
of  the  magnitudes  of  the  potential  propagation  effects  of  turbulence. 

In  Section  2  we  develop  models  for  the  size,  temperature,  and  velocit; 
of  the  rising  plume  as  a  function  of  altitude.  In  section  5  we  relate 
the  plume  parameters  to  the  strength  of  the  plume  turbulence  •  a  I 
develop  models  for  optical  signal  propagation  through  this  t;n '"•jIot- ; 
Sample  calculations  and  results  are  presented.  Th:  mode’s  dec e 1 ?pe 
here  can  be  used  in  designing  experimental  plume  tests  In  tarn,  ; h  ■ 
future  experimental  test  results  will  provide  the  needed  data  :,;.se  tV; 
verification  and  improvement  of  the  present  models.  The  plume  modi j 
development  will  be  an  iterative  process  between  the  theoretical  .;;j 
experimental  phases. 


SECTION  2 
PLUME  MODELS 


The  turbulent  buoyant  plume  has  been  an  area  of  research  for  more 
than  40  years  (References  1  to  6).  Various  techniques  have  been  used. 
Zel'dovich,  Rouse,  and  Batchelor  used  a  similarity  approach.  Schmidt 
and  Taylor  made  assumptions  about  the  basic  turbulent  flow  processes. 
Morton,  Taylor,  and  Turner  were  able  to  analyze  nonsimilar  flows  by  as 
suming  the  ambient  air  entrainment  rate  is  proportional  to  the  local 
mean  velocity.  Normally,  the  loss  of  energy  from  the  plume  by  radia¬ 
tion  is  ignored.  In  Reference  7  Morton  examined  the  radiation  loss 
and  derived  approximate  formulas  for  its  effect. 

Figure  1  shows  a  sketch  of  an  ideal  fire  plume.  We  will  solve  for 
the  plume  radius,  temperature,  and  rise  velocity  as  a  function  of  alti 
tude.  In  our  analysis  we  will  make  the  following  assumptions: 

•  The  plume  is  in  pressure  equilibrium  with  the  ambient 
atmosphere  at  all  altitudes  along  its  length. 

•  At  a  given  altitude  the  plume  has  a  radial  "top  hat" 
profile.  That  is,  the  temperature,  density,  and 
velocity  are  constant  across  the  plume  section. 

•  The  entrained  ambient  air  is  instantaneously  mixed 
with  the  plume  air. 

•  The  radiation  loss  from  the  plume  is  negligible. 

These  are  typical  assumptions  of  plume  analysis.  Experimental  measure 
ments  indicate  that  in  well  developed  plume  flows  the  temperature  and 
velocity  tadial  distributions  are  Gaussian,  with  the  maximum  values 
occurring  on  the  plume  centerline.  But  since  all  plume  formulas  are 
in  terms  of  the  integral  of  the  profiles  over  the  radial  cross  section 
only  the  averaged  temperature  and  velocity  values  are  determined. 
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Later  we  will  relate  the  "top  hat"  profile  values  to  those  of  the 
Gaussian  profile. 

We  adopt  the  Morton  (Reference  7)  entrainment  hypothesis  for 
buoyant  plumes. 


where 

a  =  entrainment  constant  at  a  given  altitude 

_3 

p  =  mean  plume  density  at  the  given  altitude  (g  cm  ) 

p  =  ambient  air  density  at  the  given  altitude  (g  cm  ^) 
a 

aQ  =  the  well  established  entrainment  constant  for  weakly 

buoyant  plumes  (0.116  for  "top  hat"  profile). 

The  entrainment  rate  of  ambient  air  mass  per  unit  plume  surface  area  is 

a  p  v,  where  v  is  the  local  plume  rise  velocity, 
a 

For  this  first  model  we  ignore  radiation  loss  terms.  Thus  our 
analysis  corresponds  to: 

•  Clean  plumes,  that  is,  plumes  that  are  not  optically 
thick  due  to  combustion  particulates, 
t  Dirty  plumes  with  relatively  low  initial  temperatures 
or  with  high  temperatures  accompanied  by  high  rise 
velocities  (in  this  latter  case  the  cooling  due  to 
rapid  ambient  air  mass  entrainment  dominates  the 
radiation  cooling) . 

At  a  later  date  we  will  add  radiation  cooling  to  the  plume  solution. 

Define  the  following  variables: 

z  =  altitude  above  reference  plane  (cm) 


T  =  plume  temperature  at  z  (K) 

p  =  plume  density  at  z  (g  cm"^) 

-2 

P  =  plume  pressure  at  z  (dynes  cm  ) 

r  =  plume  radius  at  z  (cm) 

v  =  plume  rise  velocity  at  z  (cm  s  ) 
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g  =  acceleration  of  gravity  (980)  (cm  s  ) 

E  =  internal  energy  (ergs  g  *) 

Mq  =  molecular  weight  of  air  at  standard  conditions  (28.964)  (g) 

Z  =  number  of  moles  of  air  at  temperature  T  and  pressure  P  per 
mole  of  air  at  standard  conditions 
Y  =  effective  gas  gamma  (ratio  of  specific  heats)  at  temperature 
T  and  pressure  P 

H  =  ambient  pressure  scale  height  at  altitude  z  (cm) 

y  7-1 

R  =  gas  constant  (8.317  x  10  )  (ergs  K  ) 

T  ,P  ,Z  ,y  are  the  corresponding  quantities  for  the  ambient  air. 

cl  cl  ci  cl 

We  take  the  equation  of  state  of  air  and  the  internal  energy  to  be 
d  ZpRT  ,  -2 

P  =  —  dynes  cm  (_) 

o 

E  =  FfiTHT  ergs  g_1  •  (3) 

The  values  for  Z  and  y  are  given  in  Tables  1  and  2  (References  8  and  9). 
For  fire  plumes,  temperatures  are  generally  below  2000  K  and  density 
ratios  (ratio  of  plume  density  to  ambient  sea  level  air  density)  are 
greater  than  0.1.  Thus  for  plumes  we  can  take 

Z  =  1 

and  use  the  following  theoretical  expression  for  y  based  on  the  quantiza¬ 
tion  of  the  vibrational  states  of  diatomic  molecules 

.  3.5  ♦  S 
Y  '  2.5  +  S 


where 


S  =  °-21  w  +  0-79  u 

w  u  . 

e  -  1  e  -  1 


2.239  x  10” 


3.352  x  io* 


The  first  term  in  the  equation  for  S  is  the  effect  due  to  oxygen  and 
the  second  term  is  due  to  nitrogen.  For  plume  temperatures  and  densities 
the  fit  to  y  is  good  to  three  significant  figures. 


14 


r^i 

>) 


o 


3 

cr 

LU 


OJ 


-Q 

*3 


p 

o 

G 

c 

A 

9 

© 

9 

9 

9 

« 

cr 

wc 

© 

© 

o 

© 

© 

© 

© 

© 

©  © 

©  © 

O 

©  © 

© 

G 

© 

© 

G 

© 

G 

O 

o 

C 

© 

O 

© 

G 

o 

© 

© 

© 

G 

© 

g* 

a 

2 

- 

1 

c 

G 

e 

A 

9 

x 

© 

9 

X 

cr 

P* 

9 

© 

9 

9 

rv 

9 

9 

•pc  cr 

cr  9 

A 

©  9 

P. 

A 

9 

© 

P“ 

A  9 

9 

9 

A 

0 

G 

A 

G 

© 

G 

e 

e 

« 

o 

O 

© 

o 

9 

X 

p 

9 

© 

cr 

p. 

p» 

p- 

© 

K 

© 

9 

—  cr 

PC  9 

X 

p.  — 

a 

Pi 

© 

p* 

9 

mm 

X 

A 

9 

© 

9 

A 

0 

© 

p 

© 

9 

X 

p 

9 

p 

p 

« 

o 

e 

o 

© 

© 

— 

a 

X 

9 

o 

— 

9 

9 

9 

9 

9 

9 

9 

9 

©  © 

P-  A 

9 

9  9 

X 

9 

X 

9 

C 

A 

© 

P 

© 

© 

— 

X 

9 

“r' 

9 

© 

O 

0 

9 

A 

“ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

IV 

IV 

K% 

*r 

•r 

*r* 

*PC 

■pi  »r 

K* 

•p 

c  a 

«  9 

9 

9  9 

X 

X 

X 

X 

X 

X  X 

9 

9 

9 

P- 

C 

E 

X 

X 

9 

9 

9 

© 

- 

- 

A 

•6 

c 

G 

o 

9 

c 

p 

9 

0 

cr 

© 

c 

© 

© 

© 

O 

© 

c 

©  © 

©  C 

O 

©  © 

c 

© 

© 

G 

O 

o 

© 

© 

© 

© 

© 

C 

© 

c 

© 

e 

© 

C 

© 

c 

© 

e 

_ 

G 

1 

o 

C 

o 

9 

9 

X. 

© 

® 

9 

x 

© 

O 

© 

IV 

9 

p- 

•pi 

P. 

9 

A/  © 

€i  9 

* 

pc  *pc 

G 

— 

A 

9 

X 

© 

X 

a 

£ 

9 

© 

X 

© 

9 

9 

0 

A 

A 

p 

G 

G 

© 

* 

c 

G 

c 

e 

© 

9 

C 

»r 

G 

© 

X 

fp  u- 

9 

P^ 

P- 

P-  ©  X 

9  — 

cr  A  A 

9 

p 

A 

c 

X 

9 

— 

9 

X 

X 

«• 

© 

9 

C 

p 

O 

p 

X 

8 

P 

* 

c 

O 

© 

G 

O 

© 

IV  P 

9 

O 

— 

cr 

© 

9 

9 

9 

9 

9 

9 

9  9 

9  O 

o 

—  A 

PC 

X 

9 

— 

A 

'** 

9 

9 

— 

9 

9 

0 

p- 

© 

9 

“ 

9 

© 

“ 

0 

X 

9 

O 

“ 

- 

- 

- 

- 

- 

- 

- 

- 

- 

n  iv  <v  r  r 

•p* 

CP' 

•PC 

«PC 

*»C  PI 

pc  9 

9 

9  9 

9 

0 

0 

0 

X 

X 

X 

X 

X 

C 

9 

9  ip 

p- 

i-* 

X 

X 

X 

9 

9 

« 

z 

A- 

X 

o 

© 

© 

P 

© 

9 

X 

»v 

9 

© 

«v 

o 

© 

© 

© 

O 

© 

© 

O 

©  o 

©  o 

O 

©  G 

o 

c 

© 

© 

© 

© 

© 

C 

© 

e 

© 

© 

c 

© 

e 

© 

_ 

© 

G 

_ 

c 

© 

C 

1 

o 

o 

© 

G 

p 

9 

x 

9 

X 

c 

cr  p- 

9 

P- 

IV 

© 

m 

^  a- 

»<C  © 

9 

PC  X 

9 

0 

X 

X 

9 

p- 

9 

X 

A 

9 

X 

X 

9 

p- 

X 

X 

X 

X 

P 

A 

G 

c 

o 

o 

o 

G 

o 

— 

©  P 

p 

© 

rv  f- 

PV 

•r 

p- 

9 

»p>  »r 

9  9 

p-  p» 

X 

9  © 

PC 

A 

X 

p 

9 

A 

PI 

A 

o 

A 

PC 

mm 

© 

X 

A 

9 

X 

mm 

c 

A 

4 

- 

<r 

c 

c 

c 

o 

©  A  A 

p 

9 

o 

— 

— 

9 

9  P- 

© 

9 

9 

9 

9  9 

9  9 

9 

9  © 

© 

o 

— 

A 

PC 

9 

o 

X 

« 

9 

C 

A 

X 

a 

P 

9 

p 

9 

© 

A 

X 

a 

X 

9 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

~ 

»V 

IV 

m* 

KC 

•pi 

•PC 

-c 

•PC  *P> 

•r  »pc 

•* 

PC  O 

9 

0 

0 

0 

9 

9 

X 

X 

X 

X 

9 

9 

9 

9 

p 

X 

© 

X 

9 

C 

z 

9 

o 

G 

c 

c 

9 

x 

IV. 

IV  IT 

P 

o 

p. 

© 

o 

© 

O 

O 

G 

C 

o 

o  c 

©  © 

O 

O  © 

© 

© 

O 

© 

G 

© 

© 

o 

© 

o 

© 

C 

O 

c 

C 

© 

C 

C 

c 

C 

c 

G 

c 

© 

1 

o 

9 

G 

o 

O  \x . 

X 

9 

© 

r- 

rv 

K* 

9 

cT 

9 

© 

IT 

©  *^c 

--  G 

A 

a  a 

9 

0 

0 

9 

9 

X 

9 

X 

9 

9 

X 

9 

PC 

0 

P 

X 

p 

9 

c 

« 

c 

o 

© 

o 

© 

c 

9 

a 

X 

9 

P- 

X 

mm 

9 

9 

9 

X 

p- 

a  © 

•-  pp 

a 

cr  x 

9 

p* 

c 

9 

0 

9 

9 

A 

X 

© 

© 

X 

9 

c 

A 

X 

A 

9 

9 

X 

« 

o 

c 

= 

G 

o 

© 

*“ 

IV 

9 

9 

9 

c 

9 

9 

9 

— 

fPC  X 

9 

®  © 

9  9 

9 

9  9 

9 

9 

9 

9 

O 

© 

A 

1“ 

A 

X 

p. 

9 

© 

— 

*** 

p 

o 

X 

p 

X 

X 

X 

p- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

~ 

- 

- 

A! 

IV 

IV 

IV 

■r* 

fPC 

*P| 

T 

•pi 

*p|  *PC 

PC  *PC 

PC 

•pc  *pc 

PC 

PC 

PC 

PC 

0 

0 

9 

9 

X 

X 

X 

X 

9 

9 

9 

9 

p 

p 

p 

p 

© 

9 

c 

p 

o 

e 

o 

9 

p 

p 

p 

9 

o 

© 

c 

© 

© 

© 

G 

O 

o 

o  © 

©  © 

© 

©  © 

C 

o 

o 

O 

© 

© 

C 

o 

© 

o 

© 

© 

© 

C 

G 

© 

.. 

C 

o 

o 

© 

© 

© 

G 

1 

o 

o 

o 

o 

a 

mm 

9 

G 

CC 

9 

© 

mm 

G 

cr 

9 

P- 

<v 

— 

9 

p* 

cr  *p\  x  © 

mm 

cr  9 

9 

o 

X 

A 

0 

X 

o 

X 

« 

9 

9 

a 

0 

o 

P 

A 

X 

9 

X 

9 

A 

X 

© 

O 

« 

c 

o 

G 

c 

o 

O 

© 

9 

O 

9 

9 

cr 

X 

9 

© 

CT 

—  l/c  9  9  *pc  ©  Aj  v 

X 

9 

A 

PC 

9 

p- 

PC 

© 

© 

A 

p. 

X 

© 

X  X  9 

0 

9 

A 

X 

« 

G  G 

G 

o 

G 

© 

— 

IV 

p 

9 

9 

o 

— 

rv 

H»l 

9 

9 

« 

9 

— 

*pc  C 

cr  9  p* 

p-  © 

© 

X 

X 

9 

9 

9 

9 

— 

PC 

— 

o 

9 

p- 

X 

G 

— 

X 

< 

© 

9 

X 

« 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

iv  iv  *v  iv  iv 

A. 

rv  rv  ^ 

•PC  *PC 

PC  PC 

PC 

pc  P* 

PC 

PC 

P> 

P» 

P 

P 

P 

0 

a 

9 

X 

X 

X 

X 

X 

9 

9 

9 

9 

9 

P 

P 

X 

9 

X 

• 

o 

o 

c 

o 

IV 

© 

© 

c 

© 

o 

© 

O 

o 

©  O 

©  ©  © 

©  C 

Q 

© 

© 

© 

O 

O 

O 

O 

© 

© 

© 

c 

© 

© 

c 

G 

© 

C 

© 

G 

© 

© 

© 

A 

o 

o 

o 

c 

O 

9 

a; 

9 

9  P 

9 

X  LP 

9 

IV  A, 

9 

9 

P»  PC 

—  9 

9 

©  P* 

o 

© 

o 

A 

X 

9 

« 

X 

© 

PC 

© 

9 

9 

9 

X 

P 

A  9 

0 

P 

t 

o 

o 

c 

o 

o 

e 

o 

o 

9 

cr 

P*  ©  ip' 

© 

X  9 

9 

©  »*c  « 

A  9 

tr  ir  A  9 

a 

X 

P. 

X 

* 

P 

0 

9  X 

r- 

p 

9 

X 

© 

9 

X 

9 

p 

9 

« 

9 

c 

o 

© 

o 

o 

©  A  A  X  « 

9 

© 

— 

IV  IV  pp  cT 

9 

P*  9  © 

A  *T 

9 

X  9 

9 

p» 

X 

X 

X 

9 

9 

— 

PC 

9 

9 

A 

0 

9 

r» 

X 

C 

*“ 

p. 

X 

X  A 

© 

X 

O 

— 

— 

- 

— 

— 

- 

— 

- 

- 

- 

— 

- 

rv  «v  iv  rv 

IV  rv 

rv  iv  rv  »pt 

PC  Pc 

P>  PC  P>  «c 

■PC 

PC 

■* 

pc 

pc 

PC 

0 

0 

0 

0 

X  X 

XXX 

9 

9 

9 

o 

** 

X  9 

A 

__ 

c 

o 

© 

© 

p 

© 

■o 

© 

iV 

9 

© 

© 

©  o 

_ 

o 

© 

© 

©  © 

©  © 

© 

©  © 

o 

© 

© 

o 

© 

© 

© 

© 

s> 

G 

© 

C 

c 

© 

c 

© 

© 

_ 

C 

© 

O 

a 

-I 

1 

o 

o 

o 

c 

c 

© 

9 

p 

G 

9 

cr  iv  ** 

l/' 

© 

IT 

cr 

p  9  P* 

©  A 

A 

9  p* 

o 

A 

9 

X 

r*. 

X 

r-  9 

G 

tr- 

9 

© 

P 

A 

X 

>— • 

o 

« 

e 

o 

o 

o 

o 

© 

9  «L 

P 

O 

rvi 

9 

P» 

9 

V 

©  9 

9 

o  iv  cr  ip 

9 

©  © 

©  ©  X 

9  © 

©  X 

P» 

9 

o 

fo 

0 

9 

0 

9 

mm 

P  X 

© 

© 

p 

mm 

<• 

V. 

« 

c 

c 

c 

c 

© 

o 

© 

— 

a 

9 

9 

© 

© 

— 

— 

IV 

*pc 

•PC 

9 

9  P» 

©  9 

©  A  p* 

0 

0  X 

9 

9  P- 

X 

X 

9 

© 

A  X  r- 

e 

A 

0 

X 

p 

X 

9 

— 

© 

P 

O 

* 

u. 

c 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

rv  rv 

IV  <V 

iv  rv 

i\iIVI\iI\iIVIVi''P*p' 

PC 

PH 

P* 

P 

P 

PC 

pc 

0 

0 

0 

0 

X  X 

X  X  X  X- 

X 

9 

P 

- 

X 

U 

Uj 

x 

cr. 

—1 

e 

o 

© 

e 

o 

9 

© 

9  ir 

© 

C 

© 

© 

© 

O 

o 

O 

o  c 

©  © 

© 

o  o 

o 

© 

e 

© 

© 

©  © 

© 

G 

© 

C 

C 

© 

© 

© 

© 

c 

© 

o 

© 

© 

c 

c 

Uj 

o 

c 

o 

o 

© 

© 

P 

9 

V 

9 

© 

p- 

B 

p* 

9 

9 

—  9 

•pi  9 

X 

»  X 

X 

X 

X 

X 

© 

PC 

IP 

A 

X 

PC 

0 

0 

A 

© 

X 

© 

9 

A 

A 

9 

X 

i 

c 

c 

o 

c 

© 

C 

P 

9 

9 

IV 

© 

IT 

tt 

IV 

P 

IP  a 

p.  ,r 

a  p» 

PC 

PC 

A 

9 

9 

PC 

A 

9 

A 

9 

e 

X 

0 

9 

9 

c 

« 

0 

a 

c 

A 

P 

X 

« 

o 

o 

9 

c 

o 

C 

o 

rv  ir 

cT 

® 

© 

© 

© 

O 

*M 

Aj  (V 

•PC  9 

cr  9  p. 

«  9 

o 

A 

A  pc 

0 

9 

» 

X 

© 

PC 

X 

P 

© 

0 

9 

3 

m— 

■u 

X 

* 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

> 

X 

o 

* 

*■* 

*■“ 

““ 

“ 

*** 

IV 

IV 

<V  V 

IV 

IV 

fV  rv  iV  IV 

A  A 

AAA 

PC 

PC 

PC 

PC 

** 

PC 

PI 

PC 

pi 

0 

0 

9 

0 

0 

X 

X 

X  X 

X 

X 

9 

p 

X 

o 

c 

9 

o 

e 

_ 

X 

_ 

o 

9 

cr 

9 

© 

o 

© 

G 

© 

© 

o 

O 

©  © 

G  © 

© 

G  © 

o 

© 

o 

O 

© 

C 

© 

e 

© 

© 

o 

© 

C 

© 

© 

_ 

© 

c 

© 

c 

© 

_ 

G 

►- 

» 

G 

o 

c 

c 

c 

o 

9 

o 

9 

mm 

O 

9 

V 

ir 

X 

9 

X 

9 

cr  9 

«  PI 

P*  X 

9 

0 

0 

9 

A 

9 

© 

© 

A 

X 

0 

P 

p 

X 

p 

X 

A 

A 

cr 

« 

o 

9 

c 

e 

© 

© 

IV  P 

9 

9 

P* 

9 

9 

•PC 

9 

9 

»v 

9 

—  9 

A  9 

9 

•PC  — 

9  IP 

X  pc  o 

C 

© 

A 

C 

o 

9 

p 

mm 

c 

p 

p 

C 

X 

X 

9 

u 

• 

o  o 

c 

9 

G 

© 

G 

mm 

•• 

IV 

9 

9 

9 

9 

© 

© 

© 

© 

mm 

V 

IV  A 

PC  P« 

9 

X  9 

9  P- 

X 

9 

© 

e 

P 

X 

9 

p. 

© 

X 

O 

p 

X 

p 

X 

G 

p 

C 

X 

— 

u. 

c 

o 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

■ 

•  • 

•  • 

• 

•  • 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

* 

• 

• 

• 

• 

• 

• 

C 

T> 

Uj 

— 

*“ 

“ 

— 

— 

— 

— 

— 

IV  IV 

IV 

iv 

iv  rv 

A  A 

A  A 

A 

AAA 

A 

A 

A 

** 

PC 

PC 

** 

PC 

PC 

p 

0 

0 

© 

0 

0 

9 

X  X 

ll 

9 

P 

z 

— • 

•— 

O' 

c 

o 

o 

e 

c 

c 

P 

9 

© 

9 

9 

o 

C 

© 

o 

c 

e 

© 

© 

©  © 

©  © 

© 

©  © 

© 

C 

© 

C 

© 

c 

© 

© 

G 

c 

C 

c 

C 

© 

c 

c 

© 

C 

© 

G 

C 

C 

»— 

X 

o 

o 

o 

c 

© 

© 

IV 

9  cr  a 

9 

ir 

X 

9 

p» 

*PC 

cr 

X 

Aj 

X 

©  — 

*  o  cr  cr  © 

9 

A 

PC 

X 

X 

© 

9 

X 

p 

A 

9 

X 

P- 

A 

9 

X 

9 

X 

P 

1 

c 

o 

o 

e 

© 

© 

9 

p 

IV 

PJ  X 

9 

9 

9 

X 

—  cr 

«  p1 

P 

A 

9 

X 

0 

X 

X 

O 

p* 

X 

0 

O 

9 

o 

ot 

P 

p 

mm 

P 

A 

9 

jt 

• 

e  e. 

c 

o 

o 

9 

e 

© 

* 

IV 

^■c 

9 

X 

9 

9 

9 

o 

© 

© 

© 

—  — 

x  A 

A  *r  pc 

0 

0 

X 

9 

9  P“ 

9 

— 

PC 

0 

X 

9 

r* 

X 

9 

c 

A 

p 

X 

9 

8 

P 

9 

a 

- 

a 

= 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

rv  ?<i  *V  rv 

A  A 

A  A 

A 

AAA 

A  A 

A  A 

A 

A 

■" 

PC 

p> 

PC 

* 

0 

0 

9 

* 

9 

0 

X 

X 

9 

X 

Uj 

cc. 

“ 

X 

O 

c 

c 

o 

o 

c 

© 

V  «v 

© 

9 

9 

— 

© 

© 

O 

© 

© 

o 

© 

© 

©  © 

o  © 

© 

©  C 

© 

O 

© 

O  C 

O 

C 

© 

© 

© 

© 

© 

c 

© 

© 

s 

e 

© 

9 

C 

c 

© 

£ 

© 

X 

« 

o 

o 

G 

<9 

© 

O 

p 

X 

9 

X 

9 

© 

P* 

P- 

ru 

9 

p* 

P> 

9 

©  »p> 

X  X 

9 

«  © 

P- 

9 

p 

©  X 

C 

X 

V 

0 

0 

9  r- 

X 

P 

9 

p 

p 

P 

p 

P 

P 

p 

2 

« 

o 

c 

9 

9 

e 

© 

c 

9 

9 

9 

9 

G 

«r 

IV 

© 

V 

cr 

X 

G 

A 

X  P- 

9  AJ 

•r 

©  A  X 

9 

P 

X 

A 

P* 

9 

0 

M 

X 

P* 

P* 

9 

X 

0 

p 

p 

0 

9 

X 

© 

0 

r 

z 

p 

o 

9 

o 

e 

© 

© 

© 

© 

e 

IV 

9 

c 

X 

9 

9 

9 

O 

C 

©  © 

C  — 

—  A 

A 

A  PC 

O 

o 

9 

X 

e 

A 

p 

0 

X 

9  P* 

X  9 

C 

p 

X 

p 

© 

o 

p 

• 

• 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

■ 

• 

• 

• 

• 

• 

« 

• 

• 

• 

* 

• 

X 

““ 

”” 

— 

*■ 

"" 

rv  A. 

A  A 

A  A 

A 

A  A 

A 

A- 

A 

A 

A 

A 

A  A  pc 

p 

PC 

P» 

P 

p 

P 

P 

p 

0 

0 

0 

0 

X 

X 

1 

9 

o 

© 

9 

© 

e 

9  «  IT 

© 

© 

o 

©  © 

O 

© 

© 

© 

O  O 

e  © 

e 

O  © 

o 

© 

© 

O 

© 

© 

O 

© 

o 

© 

© 

© 

» 

© 

© 

o 

© 

© 

G 

© 

O 

e 

© 

c 

p 

« 

G 

e 

a 

C 

c 

© 

IV. 

IV 

IV 

P» 

9 

9 

p 

9 

© 

mm 

cr 

P- 

o  Pc 

cr 

X  X 

9 

PC 

A 

mm 

m— 

X 

9 

A 

PC 

© 

A  C 

9 

r«. 

A  X 

© 

P 

X 

9 

p 

p 

X 

« 

© 

o 

o 

9 

© 

© 

©  IV 

P 

IV 

Pfc 

-v  cr 

mm 

©  cr 

mm 

9 

«• 

cr  9  a 

lT  P-  9 

—  p 

X 

P. 

9 

mm 

PI 

X 

pc  A 

A 

— 

© 

©  X 

A  9 

9 

A 

9 

X 

A 

9 

p 

X 

A 

o 

O 

o 

«• 

G 

c 

© 

o 

© 

o 

* 

“* 

«v 

« 

«r  cr 

9 

ip 

• 

V 

«  X  9 

9  9 

9 

c  © 

© 

O 

© 

* 

— 

A 

9 

X 

9 

9  P* 

X 

X 

9 

© 

© 

*" 

A  A 

9 

3 

0 

AAAAAAAAAAAAAAAAAA***-“i-',#'*r  p  9  0 


J 

* 


f 


t 


a  c 

A  UJ 

p  ^ 

X  W 
ft  a. 


e  e  n  c  o  o 

c  o  n  o  e  o 

-  O  \  o 

—  p  p  A 


000©00©Oeo009GOOOO©0000  900000  OOGOOOOCOC  909  909 

o  o  o  o  c  oeccoococ  oc  c  c  occooooocoe  ooooceoooeooccoo 
oe-ooooox  9/0/orf'Cir  cir>ci/'O»roir  ©  x  ©©©e©©e©©©oe©ee©©©© 

^•l/'«'*eoC  -  —  AAPP  00XX99PPX  1  ff  a  O  <V  O  9  ©  C  A  a  o«  c  \  »  «  <  c  o  o  o 

9  9  9  9  iT  9  p  * 


15 


We  first  develop  the  plume  equations  as  a  function  of  time  and  then 

as  a  function  of  altitude.  At  some  initial  time  tQ,  consider  a  cross 

section  of  the  plume  at  altitude  zq,  with  plume  radius  rQ,  rise  velocity 

v  ,  and  of  thickness  6z  .  At  a  later  time  this  volume  element  will  be 
o  o 

at  altitude  z,  have  radius  r,  rise  velocity  v,  and  its  thickness  will  be 

6z  =  6z  v/v  .  The  mass  in  the  volume  element  at  altitude  z  is 
o  o 

M  =  7tr2p<5z  =  ttt^pSz  —  g 
p  o  v  6 

o 

The  mass  entrainment  rate  is  given  by  the  Morton  assumption 


ir  =  otp  v  (2irrSz) 
dt  a  1/2 

fP,  ' 


2a 


B)  ? 


g  s 


-1 


where  the  quantity  2Trr6z  is  recognized  as  the  surface  area  of  the  volume 
element.  The  momentum  conservation  equation  is  simply 


d  (M  v ) 
dt 


-2 

g  cm  s 


where  the  quantity  on  the  right-hand  side  of  the  equation  is  the  buoyant 
force  on  the  volume  element.  The  temperature  equation,  accounting  for 
the  internal  energy  of  the  entrained  mass,  the  expansion,  and  the  chang¬ 
ing  state  of  the  plume  air,  is 


■  ,  Y  Z  T 

Y-l  a  a  a 

y  -1  Y  Z  T  ’ 


I  dM 
M  dt 


1 


_  d(Y-l)  _  I  dZi 

Y(Y-l)  dt  Z  dt 


K  s 


-1 


At  plume  temperatures  the  dZ/dt  term  is  negligible  and  the  d(Y-l)/dt 
term  is  minor. 


Naturally,  the  properties  of  the  plume,  T,  p,  r,  and  v  do  not  depend 
on  the  volume  element  considered,  but  are  only  functions  of  altitude. 
Define  the  scaled  mass  parameter 


m 


2 

pr  v 


g  s 


(4) 


then  the  equations  can  be  rewritten  in  terms  of  m,  without  reference  to 
a  volume  element.  Thus  in  the  time  domain  the  set  of  equations  to  be 
solved  are 


i/: 


dm 

l\  mv 

dt 

I  r 

dv 

dt 

V  2  /P  V 

.)  g  -  2a  — (— ) 
/.  o  r  \  p  / 

dT 

Jpm 

\  1  dP  Ty-i 

dt  " 

Tw 

/  P  dt  |_ya-l 

1/2 


2  T 
a  a 

2  T 


1_  dm 
m  dt 


1  d(Y-l) 
Y(Y-l)  dt 


1  dZ  I 

2  dt  j 


with 


dz 

dt 


=  v 


Z  T 
a  a 

p  =  pa  Z  T 


I  d£ 

P  dt 


v 

H 


The  last  auxiliary  equation  comes  from  assuming  that  over  the  altitude 
range  of  the  plume  the  ambient  pressure  is  given  by 

-2 


(Z-Zo)/H 

P  =  P  e  v  dynes  cm 

a  ao 


and  by  the  pressure  equilibrium  assumption 
P  =  P 

a 

In  an  established  plume  with  constant  heat  injection  rate  the  plume 
parameters  are  only  a  function  of  altitude,  independent  oi  time.  We  re¬ 
write  the  equations  eliminating  the  time  variable  by  the  substitution 

dt  =  ^  . 

v 

The  set  of  plume  equations  then  becomes 

1/2 


dm  _ 
t—  =  2a 
dz  o 


ft)  T 

dv  .  (^a  _  jW  (fa\ 

dz  \  p  /v  r  V  P  / 


(5) 


(6) 
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£T 

dz 


Y-l 


L.V1 


a 

Y 


Z 

a 

Z 


T 

a 

T 


1  dm 
m  dz 


+  1  d(Y-l)  _  I  dZ 

Y(Y-l)  dz  Z  dz 


with 


i  dP 

P  dz 


P  * 


1_ 

i 

P 

Z  T 


_a  _a 
Z  T 


(7) 


(8) 

O) 


Given  the  initial  plume  parameters  Tq  (or  equivalently  pq) ,  tq  and  vq 
at  zq,  the  equations  can  be  solved  for  the  plume  parameters  at  all  alti¬ 
tudes.  The  initial  radius  of  the  plume  is  determined  by  the  size  of  the 
fire.  The  initial  temperature  and  velocity  are  related  to  the  burning 
rate  (energy  release)  of  the  fire.  Consider  the  following  ideal  process. 
Let  air  be  at  temperature  Tj,  have  density  and  be  at  pressure  P.  Add 
heat  of  amount  q  (ergs/g)  and  let  the  air  expand  at  constant  pressure. 

The  energy  equation  for  this  process  is 


E2  = 


Ej  ♦  q  - 


dV 


-1 

ergs  g 


Using  Equations  2  and  3  and  solving  for 

v1  qM<. 


Z1  Y1  V1 

2  Z2  Y2  Yj-1  *1 


Z2R 


Given  T^,  Y^»  q»  and  P,  Equation  10 

Y2  by  iteration  using  the  y,  Z  tables, 
sumptions  that  at  plume  temperatures 


T2  we  have 

K  .  (10) 

can  be  solved  for  T2,  Z2>  and 
If  we  make  the  simplifying  as- 


Yi  s  Y2  =  Ya  =  1-4  for  air 

then  the  equation  reduces  to  the  usual  constant  pressure  equation 

i 

i 

i 
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We  have  written  a  computer  routine  which  solves  the  set  of  plume 
Equations  S  through  9.  The  inputs  to  the  routine  are  the  initial  plume 
parameters  Tq>  tq,  vq  and  the  ambient  atmospheric  profile.  The 
atmospheric  profile  is  specified  by  giving  the  ambient  air  density  at 
the  ground  and  the  ambient  temperature  as  a  function  of  altitude 
(either  as  a  set  of  tabular  values  or  as  a  constant  lapse  rate) .  The 
ambient  air  density  and  pressure  at  any  point  can  then  be  calculated 
assuming  hydrostatic  equilibrium. 


The  outputs  of  the  plume  routine  are  the  "top  hat"  profile  values 
of  plume  radius,  temperature,  velocity,  and  density  as  a  function  of 
altitude.  These  "top  hat"  values  represent  the  radially  averaged 
values  of  the  actual  profile.  The  initial  (just  above  the  flame  tips) 
radial  profiles  are  determined  by  the  fire  configuration.  As  the 
plume  rises  and  mixes  with  ambient  air,  the  plume  temperature  and 
velocity  rapidly  approach  Gaussian  distributions.  We  now  consider 
the  relationship  between  the  "top  hat"  and  Gaussian  profiles. 


Let  the  "top  hat"  parameters  calculated  by  the  computer  program  at 
a  given  altitude  be 

T  =  plume  temperature  (K) 
v  =  plume  rise  velocity  (m  s"1) 

—  _3 

p  =  plume  density  (g  cm  ) 

R  =  plume  radius  (m)  . 


Take  the  radial  distribution  of  the  temperature  and  velocity  as 


(14) 


(15) 


where 

r  =  radial  distance  from  the  plume  centerline  (m) 

T(r)  *  plume  temperature  at  r  (K) 
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Ta  =  ambient  temperature  (K) 

Tc  =  plume  centerline  temperature  (K) 

o  =  standard  deviation  of  radial  temperature  distribution  (it 

1  -1 
v(r)  =  plume  rise  velocity  at  r  (ms  ) 

V£  =  plume  centerline  rise  velocity  (m  s  *) 

av  =  standard  deviation  of  radial  velocity  distribution  (m) 

Note  that  it  is  the  excess  temperature  (temperature  above  ambient) 
that  is  Gaussian  distributed.  Using  the  assumption  of  pressure  equi¬ 
librium,  the  density  relationship  is  simply 

p(r)T(r)  =  pT  =  p  T  (16) 

a  a 

where 

p(r)  =  plume  density  at  r  (g  cm  3) 
p  =  ambient  air  density  (g  cm  3)  . 

3. 

The  experimental  measurements  of  George,  Alpert  and  Tamanini, 
Reference  10,  indicate  that  the  relation  between  the  "top  hat"  radius 
R  and  the  Gaussian  radial  standard  deviation  in  temperature  is 


The  measurements  of  Reference  10  found  that  the  velocity  profile  was 
wider  than  the  temperature  profile  by  about  10  percent.  The  best 
previous  experimental  measurements  were  those  of  Rouse,  et  al,  Ref¬ 
erence  4,  which  found  the  velocity  profile  to  be  narrower  than  the 
temperature  profile  by  30  percent.  The  conclusion  in  Reference  10 
after  analysis  of  the  various  experimental  measurement  techniques 
is  that  to  within  experimental  error  the  velocity  and  temperature 
profiles  can  be  assumed  to  be  equal  in  width.  Thus  we  set 


a 

v 


m 


(18) 
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In  terms  of  R  ,  the  profiles  become 

(19) 

(20) 

For  the  "top  hat"  parameters  to  correspond  to  the  averaged  Gaus¬ 
sian  parameters  requires  that  the  mass,  momentum  and  buoyancy  of  the 
plume  cross  section  be  the  same  for  both  profiles. 


(21) 

(22) 


a)  dr  <23> 

The  buoyancy  Equation  23  reduces  to  the  mass  Equation  21  and  thus  con¬ 
tains  no  additional  information.  Substituting  Equations  16  and  19 
into  the  mass  Equation  21  and  integrating  we  have 


buoyancy: 


ffR2(P  -  Pa)  = 


/R 

27rr(p(r)  - 


(24) 


(25) 
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The  excess  centerline  temperature,  ,  is  a  factor  of  2.31  times 

the  excess  "top  hat"  average  temperature,  T  -  T  ,  at  low  excess  tern- 

Si 

peratures  and  rises  to  a  factor  of  2.97  when  T  -  reaches  1000  K . 

Substituting  Equations  16  and  20  into  the  momentum  Equation  22 
and  solving  we  have 


Hence  the  ratio  of  to  v  is  the  same  as  the  excess  temperature 

ratio,  ranging  from  2.31  to  2.97  as  T  -  T  ranges  from  0  to  1000  K. 

We  have  used  the  plume  computer  routine  to  solve  for  the  "top  hat" 

profile  parameters  for  a  variety  of  initial  conditions.  Then  assuming 

that  the  actual  profile  is  Gaussian  at  all  altitudes,  we  have  used 

Equations  25  and  26  to  convert  from  "top  hat"  to  Gaussian  parameters. 

Figures  2  through  28 show  the  results.  Plotted  are  the  plume  center- 

line  excess  temperature,  T  -  T  ,  the  centerline  velocity,  v ^  , 

c  a.  c 

and  the  "top  hat"  radius  R  as  a  function  of  altitude.  The  radius 
R  represents  the  2-sigma  point  of  the  Gaussian  temperature  and  veloc¬ 
ity  distribution.  Thus  the  temperature  and  velocity  decreases  from 
the  centerline  values  at  r  =  0  to  values  a  factor  of  e  =  0.135 
of  the  centerline  values  at  r  =  R  .  The  altitude  in  the  figures  is 
the  altitude  above  the  reference  plane  where  the  initial  plume  param¬ 
eters  are  specified;  it  can  be  thought  of  as  the  altitude  above  the 
tips  of  the  flames.  The  ambient  atmosphere  was  assumed  to  have  a 
ground  level  temperature  of  288  K,  and  have  a  constant  neutral  stability 
lapse  rate  (dT  /dz  =  -0.0098  Km’).  Until  the  plume  temperature  ap- 

a 

proaches  ambient  temperature,  the  plume  parameters  are  relatively  in¬ 
sensitive  to  the  lapse  rate.  For  the  short  range  of  altitudes  shown 
in  the  figures,  the  plume  parameters  would  be  essentially  the  same  for 
any  reasonable  lapse  rate.  As  the  plume  approaches  ambient  conditions, 
the  plume  parameters  become  highly  dependent  upon  the  temperature  pro¬ 
file  (lapse  rate). 
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The  initial  plume  radii  for  the  sample  cases  are  taken  as  1,  5,  and 
10  meters.  These  are  chosen  to  cover  the  likely  range  of  individual 
battlefield  fires.  A  1-m  fire  corresponds  approximately  to  a  piece  of 
equipment  or  small  vehicle  such  as  a  jeep;  a  5-m  fire  represents  large 
vehicles  such  as  trucks  or  tanks;  and  a  10-m  fire  is  representative  of 
structures  such  as  houses.  Area  fires  such  as  burning  fields  cannot  be 
represented  by  a  single  fire  plume.  The  initial  plume  centerline  tem¬ 
peratures  are  taken  as  500,  1000,  and  1500  K  above  ambient  temperature. 
These  temperatures  correspond  to  cool,  normal,  and  hot  fires.  For  each 
initial  radius  and  temperature  combination,  three  initial  centerline 
velocities  are  chosen,  corresponding  to  strong,  medium,  and  weak  ac¬ 
celeration  of  the  heated  air  at  the  base  of  the  plume. 

We  have  specified  the  initial  ambient  and  plume  temperatures  rather 
than  the  initial  air  densities.  By  specifying  the  temperatures,  the 
results  for  the  mean  plume  parameters  are  essentially  independent  of 
the  elevation  of  the  ground  surface.  The  mean  results  apply  equally 
well  to  sea  level  as  to,  say,  the  4000-foot  elevation  at  the  White 
Sands  Missile  Range.  The  independence  of  ground  elevation  (initial 
ambient  density)  is  readily  apparent  upon  examination  of  the  plume 
equations  5  through  9.  At  plume  temperatures  and  densities 

z  =  1 


The  plume  density  enters  only  weakly  in  the  determination  of  y ;  and  at 
plume  temperatures  and  densities,  y  differs  only  slightly  from  the  ambient 
value  of  1.4. 


25 


ALTITUDE  (m) 


VELOCITY  (ms'1) 


Figure  3.  Plume  parameters  for  initial  conditions  T  -T  ■ 500  K.  V  =5  m  s'1. 
R  *  1  m.  c  a  c 
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Figure  5.  Plume  parameters  for  Initial  conditions  Tc-T  *  1000  K,  Vc *  2  m  s 
R  =  1  m. 
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Figure  6.  Plume  parameters  for  initial  conditions  T  -T  =1000  K,  V 
R  =  1  m.  c  a  c 
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Figure  7.  Plume  parameters  for  Initial  conditions  T  -T  *  1000  K,  V  =  10  m  s"1, 
R«  1  m.  c  a  c 
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Figure  11.  Plume  parameters  for  initial  conditions  T.-T  *500  K,  V  =5  m  s”1, 
R*5  m.  c  a  c 
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Figure  12.  Plume  parameters  for  initial  conditions  T  -T  =500  K,  V  =10  m  s"1 
R=5  m.  c  a  c  ’ 
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Figure  13.  Plume  parameters  for  Initial  conditions  T  -T  *500  K,  V  =15  m  s’1, 
R*5  m.  c  a  c 
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Figure  15.  Plume  parameters  for  initial  conditions  T  -T  =1000  K,  V  =10  ms1, 
R=5  m. 


39 


ALTITUOt  (ni) 


0 


10 


VELOCITY  (m  s*1) 


20 


30 


Figure  16.  Plume  parameters  for  initial  conditions  T  -T  =1000  K,  V  =15  m  s'1, 
R=5  m.  c  a  c 
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Figure  19.  Plume  parameters  for  Initial  conditions  T  -T  =1500  K,  V  =15  m  s"1, 
R=5  m.  c  a  c 
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Figure  20.  Plume  parameters  for  Initial  conditions  T  -T  =500  K,  V  =10  r.i  s"1, 
R*10  m.  c  a  c 
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Figure  21.  Plume  parameters  for  Initial  conditions  T  -T  =500  K,  V  =15  m  s'1. 
R-10  m.  c  a  c 
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Figure  22.  Plume  parameters  for  Initial  conditions  T  -T  «500  K,  V  =20  m  s"1, 
R-10  m.  c  a  c 
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Figure  23.  Plume  parameters  for  Initial  conditions  T  -T  =1000  K,  V  *10  m  s'1. 
R*10  m.  c  a  c 
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Figure  24.  Plume  parameters  for  initial  conditions  T  -T  *1000  K,  V 
R*10  m.  c  a 
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Figure  25.  Plume  parameters  for  initial  conditions  T  -T  =1000  K  V 
R*10  m.  c  a 
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Figure  26.  Plume  parameters  for  initial  conditions  T  -T  *1500  K,  V  *10  m  s'1 
k-iu  m.  c  a  c  * 
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Figure  27.  Pluroperemeters  for  initial  conditions  Tc-Ta.1500  K,  V£=15  *  s'1. 
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Figure  28.  Plume  parameters  for  Initial  conditions  T  -T  =1500  K,  V  =20  m  s'1 
R=10  m-  c  a  ’  c  » 


SECTION  3 

PROPAGATION  MODELS 


INTRODUCTION 

The  index  of  refraction  gradients  and  fluctuations  within  a  tur¬ 
bulent  hot  fire  plume  can  produce  significant  effects  upon  the  propaga¬ 
tion  of  electro-optical  signals.  There  are  two  classes  of  effects; 
the  first  is  due  to  the  mean  properties  of  the  plume  and  the  second 
is  due  to  the  turbulent  fluctuations  about  the  mean.  The  radial  and 
vertical  mean  temperature  gradients  in  the  plume  will  produce  corre¬ 
sponding  mean  index  of  refraction  gradients.  These  mean  index  gradients 
will  cause  an  electro-optical  ray  path  to  bend  away  from  the  straight 
line  propagation  path  towards  the  regions  of  higher  index.  The  turbu¬ 
lent  fluctuations  of  the  index  of  refraction  about  the  mean  index  will 
cause  random  propagation  effects.  These  turbulent  effects  include  beam 
spreading  and  jitter,  signal  decorrelation,  and  phase  and  amplitude 
scintillations.  First  we  will  examine  the  mean  bending  due  to  the 
thermal  gradients  and  then  the  turbulent  propagation  effects. 

It  is  the  spatial  and  temporal  distribution  of  the  index  of  refrac¬ 
tion  which  determine  the  electro-optical  propagation  effects.  Clean 
air  is  essentially  nonabsorbing  at  visible  optical  wavelengths  and  for 
microwaves  at  frequencies  less  than  about  72  GHz.  The  index  of  refrac¬ 
tion  is  then  real  and  is  given  by  (References  11  and  12} 

.-1  ■  (l  ♦  «  J°~3).  »*],  (27) 

where 
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P  =  atmospheric  pressure  measure  in  millibars 
(1  standard  atmosphere  =  1013.2  millibars: 

A  =  radiation  wavelength  (pm) 

e  =  water  vapor  pressure  (millibar?) 

T  =  absolute  temperature  (K) 

(-0.120  at  optical  frequencies 

B’ 

14810  at  microwave  frequencies. 

At  optical  frequencies,  the  humidity  dependent  term  is  nog 1 . r  ie . 

In  the  infrared  and  for  the  higher  microwave  frequence"? ,  ibsc 
tion  effects  enter  and  anomalous  dispersion  i  varving  dependence 
wavelength)  occurs.  Then  the  index  of  refraction  can  be  a  raridi. 
changing  function  of  frequency  and  the  simple,  nondispei si ve  ir.oex 
of  refraction  formula  no  longer  applies.  Normally,  however,  trie  nor  - 
dispersive  formulation  is  adequate  in  the  normal  atmospheric  iiaru. - 
mission  windows,  where  systems  are  designed  to  operate. 

MEAN  BENDING  EFFECTS 

Consider  the  propagation  geometry  shown  in  Figure  29.  A  horizonta 
ray  path  passes  through  a  fire  plume  at  altitude  h  at  a  distance  d 
from  the  plume  centerline.  The  cross  section  of  the  plume  at  altitude 
h  is  shown  in  Figure  30.  For  calculation  purposes  we  place  a  carte¬ 
sian  coordinate  system  with  origin  on  the  plume  centerline  at  altitude 
h  ,  with  the  z  axis  vertical  (coincident  with  the  plume  centerline), 
x  axis  horizontal  and  parallel  to  the  ray  path,  y  axis  horizontal  and 
perpendicular  to  the  ray  nath. 

The  mean  radial  and  vertical  index  of  refraction  gradients  in  the 
plume  will  cause  the  ray  path  to  bend  away  from  the  geometrically 
straight  ray  path.  The  radial  gradient  will  cause  a  bending  in  the 
y  direction  away  from  the  plume  centerline  and  the  vertical  gradient 
will  cause  a  vertical  bending  in  the  positive  z  (upward)  direction. 

The  total  mean  angular  deviation  caused  by  a  single  plume  is  expected 
to  be  small,  on  the  order  of  a  milliradian.  But  the  total  linear 


Figure  29  .  Ray  path  geometry 


Figure  30.  Cross  section  of  ray  path  through  pi 


offset  from  the  straight  line  path  at  the  receiver  depends  on  the 
distance  between  the  plume  and  the  receiver.  At  1000  meters  beyond 
the  plume,  a  one-mi  1 1 iradian  angular  bending  will  produce  a  linear 
offset  of  one  meter.  Whether  this  magnitude  of  angular  bending  is 
significant  depends  upon  the  particular  electro-optical  system. 


Define 


U 


dx 

ds 


w 


dz 

ds 


where  s  is  the  distance  along  the  true  ray  path.  "h>n  ri\ 
for  rav  tracing  are 


dU 

ds 


j. 

n 


9n 

3x 


dV 

ds 


_1 

n 


(1 


r 


(2: 


(.  30  > 


In  the  region  of  the  plume,  the  deviation  of  the  true  ray  path  from 
the  straight  line  path  will  be  extremely  small  so  we  can  set 


s  =  x 
U  =  1 

and  drop  Equation  28  from  consideration.  In  the  nondispersive  wave¬ 
length  regime, 

n  -  1  <  <  1 

so  we  can  take 
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-  =  1  . 
n 

Since  the  bending  is  so  slight 

and  we  neglect  all  second  and  third  order  products  in  Equations  29 
and  30.  Our  bending  equations  become  simply 

dV  =  3n 
dx  3y 

3W  _  3n 
3x  ~  3z 


For  x  beyond  the  plume 

V(x)  =  ?  =  tan  0  a:  0 
dx  y  y 


since  8  <<  1 

x 


W(x)  =  ~  «  tan  0  2=  0 
dx  z  z 


where  0  and  0  are  the  total  angular  bending  in  the  horizontal 

y  z 

and  vertical  planes  due  to  the  plume  mean  thermal  gradients.  Hence 


dx  rad 


r  a 

J  3z 


where  xR  ,  x^,  are  the  receiver  and  transmitter  x  coordinates. 
Since  we  are  integrating  over  a  Gaussian  temperature  distribution 
and  the  integration  limits  are  well  beyond  the  plume,  the  finite 
limits  can  be  replaced  by  ± ,  which  will  simplify  the  evaluation 
of  the  integrals. 
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Differentiating  Equation  27  for  the  index  of  refraction  we  have 


3n  _  9 (,n-l) 
9y  9y 


9n  =  3 (n- 1 J 
9z  3z 


10 


9T 

3y 


to 


91 

3t 


(33) 


Assuming  optical  and  infrared  frequencies,  the  humidity  dependent  tern 
has  been  neglected.  The  plume  is  assumed  to  be  in  pressure  e-uni  lib¬ 
rium  so  there  are  no  pressure  gradients  in  the  horizontal  •.!  .it  *.;■•!. 
There  is  a  vertical  pressure  gradient,  9P/9z  ,  the  hydros! j-  •  pi 
sure  lapse  rate;  but  the  pressure  gradient  term  is  less  than 
of  the  temperature  gradient  term  and  has  been  neglected 

The  plume  Gaussian  temperature  distribution.  Equation  19,  cur, 
written  as 


•>  .x 


2  2 
y 


T  ( x ,  y ,  z )  =  T  (z)  ♦  [T  (z)  -  T  (z)Je  ’  K 

a  C  a 


(33 


where  we  have  explicitly  indicated  the  coordinate  dependencies.  Along 
the  integration  path  y  =  d  .  Thus 


2  ,2 


3T(x,d,t)  = 
9y 


4d 


R2(z) 


9T(x,d,z)  (9Tc(z) 

- +  -Ta(*)] 


9z 


I 


7  X 

7 

-  T  (z)] 

e 

a 

r  2  ,21 

x  +  d 

9R(z) 

U3mJ 

3z 

1 36) 


K  m 
2 

x“  +  d*~ 

R2{Z)  k  m'1 
K  m 


(37) 


[n  Equation  37  we  have  neglected  the  ambient  temperature  lapse  rate, 
3T  (z)/3z  ,  in  comparison  with  the  other  terras. 

3 


The  integrals  defining  0  and 


0^  are  obtained  by  combining 


Equations  31  through  37.  These  are  non-integrable  analytically,  ex¬ 
cept  in  limiting  cases,  and  must  be  solved  numerically  in  the  general 
case.  We  will  rearrange  the  integrals  for  ease  in  graphical  and 
numerical  solution.  Define  the  following  parameters 


w  =  =  dimensionless  path  integration  parameter 


F  =  — - _ 

d  R(z) 


=  dimensionless  point  of  closest  approach 


FT  = 


T  (z)  -  T  (z) 


Ta(z) 


=  dimensionless  excess  temperature 
on  plume  centerline 


Pq  =  1013.2  millibars  =  standard  atmosphere 
Tq  =  288  K  =  typical  ground  level  temperature 

77.6  P  x  10-6  _4 

n  -1  =  - 2 -  *  2.73  x  10  =  typical  ambient  air  index 

0  o  of  refraction  minus  one  at  sea  level. 


Define  the  following  normalized  integrals 


I  (F.)  =  8  x  io3 

o  a 


IiCFd»FT)  =  8  *  1q3 


^^d’V  =  32  x  iq3 


Then  the  mean  bending  formulas  with  the  angles  expressed  in  milli- 
radians  become 
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)  milli radian. 


y2  3R  r  r  r  r  ,  i  §  ,  ,• 

fd  3; J  IlfFd,rT)  4  3z  :(  d' 


4T 

R  ‘  e 
4 ( T  -T  1  3; 

C  :j 


C  F  j  **  ^  ^  I  j » J  -|’l  mi  1 1  x  r cii.i  i  nn  s  ^  ^  «- 


For  most  wavelengths,  the  wavelength  dependent  term,  7.32  ■ 

can  be  neglected.  Note  that  the  bending  does  not  depend  on 

ical  size  of  the  plumes.  Large  and  small  plumes  which  are  > 

produce  the  same  bending.  The  similarity  is  obvious  in 

lation.  In  the  0  formula  the  term 
z 


T  -  T  3z 
c  a 

represents  the  fractional  change  in  excess  centerline  temperature  pt • 
change  in  height  of  one  plume  radius.  The  3R/3z  term  is  the  same 
for  similar  plumes.  In  fact,  in  the  weakly  buoyant  rise  phase  of  the 
plume  all  model  plumes  have 


-  0.18 
3  z 


This  0.18  model  value  is  also  confirmed  by  the  experimental  measure¬ 
ments  of  Reference  10. 

The  plume  radius,  centerline  temperature  and  gradients  are  deter¬ 
mined  by  the  plume  results  of  Section  2.  We  evaluate  the  normalized 
integrals  in  the  following  manner.  The  I  integral  has  the  analytic 
solution 


x  -2Fa  ’2¥a 
I  (F.)  =  8  x  io3(n  -1)  —  e  d  =  1.37  e  d 
o  d  o  m 
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The  1^  and  I ^  integrals  have  to  be  evaluated  numerically.  They 
are  well  behaved  integrals  and  any  standard  numerical  integration 
technique  can  be  used  for  evaluation.  For  user  convenience  we  present 
the  solutions  in  graphical  form.  Define  the  following  ratios 

It(F,,F  ) 

Fio(Fd’V  3  no3 

o  d  T 
(F , j  F  ) 

F  fF  F  )  =  — _ = _ L- 

^21LW  I1(Fd,FT)  . 

Figures  31  and  32  show  the  ratios  as  a  function  of  F,  and  F^  . 

d  T 

Naturally 

*2Fd 

VW  =  VW  Vd’V  =  l'i7  e  F10^Fd'FT)  (44) 


VW  3  VW  Wt} 


(45) 


As  an  example  calculation  consider  a  plume  with  the  following  parameters 
\  =  1  pm 

P  =  P 

o 

T  =  T 

o 

Fd  =  0-5 

F  =  WOO  =  3  47 
T  T 


R  3Tc 


T  -T  3z 
c  a 


=  -0.5 


'  I* 

3z 


-0.1  (ie,  we  are  in  the  strongly  buoyant  accelerating 
rise  phase). 


From  Figures  31  and  32  we  find 

F,  JO. 5,  3.4-)  =  0.21 
10 

Fm(0.5,  3.47)  =  2.20  . 

Using  Equations  44  and  45 

I  JO. 5,  3.47)  =  0.17 
Io(0.5,  3.47)  =  0.38  . 

Using  Equations  41  and  42 

ft  =0.30  milli rad 

y 

h  =0.12  mil li rad 

z 

TURBULENT  PROPAGATION  EFFECTS 

The  turbulence  in  a  hot  fire  plume  causes  intense  random  index  - 
refraction  fluctuations  which  can  cause  significant  effects  on  the 
propagation  of  electro-optical  signals.  These  effects  inc  ludv  beam 
spreading  and  jitter,  loss  of  coherence,  and  phase  and  amplitude 
scintillations.  We  first  consider  the  index  of  refraction  fluctua¬ 
tions  and  then  the  propagation  effects. 

Define 


A  =  77.6P  1  ♦ 


7,52  x  10 


Then  at  optical  frequencies 


,  A 
n-1  =  f 


and  the  variation  in  the  index  of  refraction 


dn  =  d(n-l)  =  A 


dP  aI l\i 
PT  +  d\T  )  | 


For  infinitesimal  variations 


-(*)■- s  • 
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At  a  point  in  either  the  natural  atmosphere  or  a  fire  plume  the  pres¬ 
sure  fluctuation  term  is  negligible  compared  to  the  temperature  fluctua¬ 
tion  term.  The  variance  of  the  index  of  refraction  fluctuations  is 
thus 


<  -  °1/T 


(49) 


where  o^T  is  the  variance  in  1/T  fluctuations.  For  weak  (ie,  small 

amplitude)  fluctuations 

2 


2  ~  T  -2 

Vr  =  7  K 


(50) 


where  a^,  is  the  variance  of  the  temperature  fluctuations. 


First  consider  a  point  in  the  natural  atmosphere.  Turbulence  will 
cause  air  parcels  originally  from  both  higher  and  lower  altitudes  to 
pass  by  the  point.  The  temperature  of  an  air  parcel  at  the  point  de¬ 
pends  on  the  mean  temperature  of  the  air  at  the  altitude  the  parcel 
originated  and  the  change  in  temperature  due  to  adiabatic  expansion  or 
contraction  in  moving  from  one  altitude  to  another.  Near  sea  level 

_3 

the  adiabatic  lapse  rate  is  -9.8  x  10  K  per  meter.  The  temperature 
fluctuations  are  small  so  that 


0 


2 


n 


2 

t4  °T 


(51) 


For  a  sea  level  pressure  of  1013.2  millibars  and  ambient  temperature 
of  288  K  at  optical  frequencies 


a2  =  9  x  10'13  o2 
n  T 


In  the  atmosphere  the  temperature  variance  associated  with  strong 

2 

turbulence  can  be  of  the  order  of  a  (degree)  so  that  the  index  of 

-12 

refraction  variance  can  be  of  the  order  of  10 

Next  consider  a  point  in  the  fire  plume.  The  temperature  fluctua¬ 
tions  in  the  hot  plume  can  be  several  orders  of  magnitude  larger  than 
the  ambient  temperature  fluctuations.  Hence  the  turbulent  index  of 
refraction  variances  can  be  four  orders  of  magnitude  larger  than  the 
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ambient  variances.  Let  T  be  the  mean  temperature  of  the  plume  at  a 
2 

point,  and  <Xj,  the  variance  ot  the  temperature  fluctuations  about  the 
mean .  Then  for 


where  T  is  the  ambient 
a 

approximation 


a" 

n 


A2a^ 

T  4 


temperature , 


we  can 


the  small  t in. tout 


But  for  larger  temperature  variations  we  cannot  ■  ‘he  »  \ 
but  must  use  the  full  variance  equation 


a 


2 


n 


ai/T 


The  temperature  fluctuation  probability  distribution  i «  renin  red 
in  order  to  solve  the  full  variance  equation.  In  well  developed  wet k 
buoyant  plume  flows  the  temperature  f luctuat.ionc  are  approximately 
Gaussian  distributed.  In  the  general  plume  flow  the  probability  dis¬ 
tributions  are  not  known.  We  will  consider  the  Gaussian  and  several 
limiting  general  flow  distributions.  We  first  consider  the  Gaussian. 
Assume  the  temperature  at  a  point  in  the  plume  has  the  probability 
distribution 


where 

a  =  normalization  constant  which  ensures  that 
2  T  -  T 

a 

f(T)dT  =  1 

i 

=  truncated  Gaussian  parameter  related  to  the  temperature 
standard  deviation  (K) . 

The  temperature  distribution  chosen  is  a  truncated  Gaussian.  The  lower 
limit  of  the  temperature  range  is  taken  as  the  natural  limit  of  the 
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ambient  temperature.  The  upper  limit  of  the  temperature  range  was 
chosen  so  as  to  have  a  symmetrical  probability  distribution  about  the 
mean.  The  upper  limit  is  the  same  number  of  standard  deviations  above 
the  mean  as  the  lower  limit  is  below. 

The  cumulative  distribution  of  the  normal  random  probability  function 
is  defined  as 


Values  of  Cq(x)  are  available  in  mathematical  handbooks;  we  also  have 
a  computer  routine  for  CQ(x).  In  terms  of  Cq(x),  the  normalization 
constant  a  is 


1 


The  temperature  variance  is  by  definition 

2  2~  —  2 

=  cn  -  or 

2  T  -  T 

a  2  —  2  2 
T  f(T)dt  -  OT  K 


(56) 


Carrying  out  the  integration  we  have  the  relationship  between  the 

2 

temperature  variance  and  the  truncated  Gaussian  parameter 


(57) 


The  defining  equations  of  a 


2 

1/T 

2 


K 


are 

-2 


(58) 
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where 


a 


Equations  59  and  oO  have  to  be  evaluated  numerical  i\  ".a-  he  •  • 

a  computer  routine  which  for  inputs  T,  "  computes  r'  ".sice 
Equations  55  through  60.  is  then  found  v-'ir-  Equation  1  - 

In  the  experimental  measurements  of  i'eit  ercc  I  >  *  »  e  r  • 

fluctuation.,  in  a  well  developed,  weakl>  vu-.nnnr  p  1  u.r,c  ♦  ;»  fci  -as 
found  that  on  the  plume  centerline 

oT  ~  0 . 4  (T  -  T  )  ■  c  ’ 

i  a. 

Figure  33  shows  the  index  of  refraction  variance  for  this  issumed 
temperature  variance.  Also  shown  is  the  variance  calculated  using  the 
small  fluctuation  approximation,  Equation  S3 .  The  simple  approximation 
equation  can  be  used  for  excess  temperatures  below  about  100  K;  at 
higher  temperatures  the  full  equations  must  be  used.  We  see  that  the 
index  of  refraction  variance  first  increases  as  the  plume  temperature 
increases,  reaches  a  peak,  and  then  declines  as  the  temperature  increases 
still  further.  The  decline  with  increasing  temperature  is  caused  by  the 
inverse  temperature  dependence  of  the  index  of  refraction.  The  peak  is 
some  three  orders  of  magnitude  above  the  typical  values  encountered  in 
strong  turbulence  in  the  atmosphere.  Note  that  in  a  plume,  the  mean 
temperature  is  a  maximum  at  the  plume  centerline  and  decreases  radially 
to  the  ambient  temperature.  Hence  even  in  a  very  high  temperature  plume 
there  will  be  a  portion  of  the  plume  at  the  temperature  of  the  maximum 
index  of  refraction  variations. 


68 


T-  Ts  (K) 


Figure  33.  Index  of  refraction  variance  as  a  function  of  excess  plume 
temperature  for  Gaussian  temperature  fluctuations. 


The  relationship  between  o^  and  mean  excess  plume  temperature  given 
in  Equation  61  is  for  well  developed,  weakly  buoyant  plume  flow.  We 
can  assume  that  the  initial  hot  plume  air  and  the  subsequently  entrained 
cool  ambient  air  are  relatively  well  mixed.  The  experimentally  measured 
Gaussian  type  temperature  distributions  at  a  point  confirm  the  well 
mixed  assumption.  But  there  are  circumstances  where  we  do  not  expect 
well  mixed  flows.  In  the  weakly  buoyant  phase  the  turbulence  is  ex¬ 
pected  to  be  of  the  dissipative  type.  In  the  low-altitude,  strongly 
buoyant  regimes  the  flow  is  primarily  convective.  Until  the  convective 
flow  decays  into  a  dissipative  flow,  mixing  may  nor.  be  well  dove  Tired . 
Also  interactions  of  the  ambient  turbulent  flow  field;  may  disrupt  '  i  .- 
smooth  mean  flow  fields  of  our  ideal  noninteracting  plume  model.  A 
likely  result  of  the  disruptions  is  the  introduction  of  volumes  of 
ambient  air  into  the  plume  flow,  or  the  breakup  of  the  plume  resulting 
in  volumes  of  the  plume  air  being  introduced  into  the  ambient  flow 
fields.  The  disruptions  will  produce  volume  elements  that  are  not  well 
mixed. 


Flows  which  are  not  well  mixed  will  have  relatively  large  tempera¬ 
ture  fluctuations  about  the  mean  temperature.  The  temperature  fluctua 
tion  probability  distribution  will  no  longer  be  of  the  Gaussian  type. 
Gaussian- like  distributions  are  characteristic  of  relatively  small 
temperature  fluctuations,  that  is,  relatively  well  mixed  flows.  For 
the  symmetric  truncated  Gaussian  the  maximum  standard  deviation  (square 
root  of  the  variance)  possible  is 


T  -  T 


°T  - 


/I 


0.  58  (T  -  TJ 

a 


T  <  T  <  2  T-T 
a  —  -  a 


Even  if  the  upper  temperature  limit  is  allowed  to  extend  to  infinity 
so  we  have  an  asymmetric  truncated  (one  side)  Gaussian  distribution, 
the  maximum  standard  deviation  is  only 


(T  -  T  ) 
a 


0.76  (T-T)  T  <  T  <  00 

cl  3 
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We  begin  our  study  of  non-well  mixed  flows  by  examining  the  limiting 
case  of  completely  nonmixed  temperature  elements. 

Consider  the  following  simple  limiting  case.  In  the  vicinity  of  a 
point  let  the  plume  air  consist  of  discrete  volumes  of  air.  The  dis¬ 
crete  volumes  are  at  either  the  ambient  temperature  T  or  some  higher 

a 

temperature  T2.  Let  f  be  the  fraction  of  the  total  volume  which  is  at 
temperature  T  .  Obviously  the  fraction  of  the  volume  at  T,  is  1  -  f  . 
The  mean  (volumetric)  temperature  for  this  simple  discrete  bimodal 
temperature  distribution  is 


T  =  faTa  +  (1  *  fa^T2  K 


and  the  variance  is 


(62) 


-  f  J2  ♦  (1  -  fft)T2  -  (T)2 

*  V^a^W2 

•  (T2-T)(T-Ta)  .  (63) 

For  discrete  volume  elements  the  index  of  refraction  variance  becomes 


o2  =  (n2)  -  [n]2 

■"af, 

A2(T2  -T)(f-Ta)  A2o2 
(TJ2)2  '  (TaT2)2 


(64) 


where 

n  =  index  of  refraction  of  air  at  temperature  T 
*  A/Ta 

n2  =  index  of  refraction  of  air  at  temperature  T2 

=  a/t2. 
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T  is  the  mean  temperature  for  discrete  volume  elements.  If  the 
volume  elements  were  uniformly  mixed  the  mean  temperature  of  the  uni 
formly  mixed  volume  would  become 


T 

u 


f  T  +  (1  -  f  )T 
am  a  am  a 


f  L  *  (1  -  f  )T 
a  2  a  a 


('>r 


where  f  is  the  fraction  of  the  total  mass  which  is  at  temperature  r 
am 

In  deriving  the  second  part  of  Equation  65,  we  used  the  pressure 
equilibrium  assumption 


p  T 
a  a 


i '  ■ 


where  p  and  p~,  are  the  densities  in  the  volume  elements  before  •  i  < . . 
a  2 

and  is  the  density  of  the  uniformly  mixed  volume. 

Next  let  us  consider  a  very  simple  model  for  the  non-well  mixed 
regions  of  the  fire  plume  and  evaluate  the  index  of  refraction  flucr  :.i 
tions.  We  assume  that  initially  at  the  base  of  the  plume  the  plume 
air  consists  entirely  of  hot  air  at  temperature  Ip.  As  the  initial 
hot  air  rises,  cold  ambient  air  of  temperature  Ta  is  entrained  into  the 
plume.  We  assume  that  a  fraction  fM  of  the  entrained  air  is  uniformly 
mixed  with  the  initial  hot  plume  air,  resulting  in  a  mixture  with 
intermediate  temperature  T,^.  The  rest  of  the  entrained  air,  of  frac¬ 
tion  (1  -  fM) ,  remains  at  the  ambient  temperature  Tg .  The  temperature 
distribution  is  thus  discrete  bimodal,  one  temperature  being  T  the 
second  Tu.  The  mixing  fraction  fu  is  a  rough  measure  of  how  well 
mixed  the  flow  is.  Figure  34  shows  a  simple  diagram  of  the  mixing 
model . 


The  volume  fractions  in  this  mixing  model  are 


f  b  =  total  volume  fraction  of  ambient  air  before  mixing 


=  volume  fraction  of  initial  hot  plume  air 

-  (1  -  f,.)f  ,  =  volume  fraction  of  air  left  at  ambient 
M  ab 

temperature  after  mixing 
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(a)  Before  mixing 


f  =  f..f  ,  =  volume  fraction  of  ambient  air  mixed  with 
ax  M  ab 

initial  plume  air. 


The  temperature  of  the  uniformly  mixed  plume  and  ambient  air  is  found 
by  applying  Equation  65  with 


f 

a 


ax 


fax  +  ^*fab) 


T 

u 


The  result  is 


TM  T„f 


T  Tr [ 1  -  f  J 
a  Fl  a  £ 


F*ax  +  V1'  fab> 


(O' 


The  mean  temperature  of  the  discrete  bimodal  volume  distribution  aftei 
mixing  is  found  by  applying  Equation  62  with 


fa  "  fa£ 


T2  =  TM 
which  gives 


T 


T  f  „ 
a  a£ 


+  V1-^ 


(68) 


Applying  Equations  63  and  64  we  have  for  the  temperature  and  index  of 
refraction  variances  after  mixing 

Cy  =  (TM-T)(T-Ta)  K2  (69) 

0  =  - - -  =  - y  (70) 

‘W  <TaTM> 


Figure  35  shows  the  index  of  refraction  variances  for  various  values 
of  the  mixing  fraction  f„.  The  initial  conditions  assumed  were  standard 


Figure  35.  Index  of  refraction  variances  for  discrete 
bimodal  temperature  distributions. 


atmospheric  pressure  and 


T  =2 88  K 
a 

T  =  1288  K 
r 

Also  shown  for  comparison  is  the  index  variance  for  the  truncated 
Gaussian  temperature  fluctuation  distribution  of  the  well  developed 
flow  with 

aT=0.4(T-Ta)  . 

As  expected,  the  variance  for  the  discrete  volumes  increases  as  the 
mixing  of  the  initial  hot  plume  air  decreases.  At  the  peak,  the  no 
mixing  case  (f^  =  0)  has  a  variance  an  order  of  magnitude  larger  than 
the  well  developed  flow  variance.  For  excess  temperatures  below  about 
100  K,  the  well  developed  flow  corresponds  fairly  closely  to  the  i‘M  - 
0.25  mixing  case.  The  sharp  falloff  as  the  excess  temperature  approaches 
1000  K  occurs  because  in  our  simple  discrete  model  initially  all  the 
plume  air  is  at  temperature  Tp  with  no  entrained  ambient  air.  Thus 
initially  the  variance  is  zero  in  the  discrete  model. 

PROPAGATION  MODELS 

We  next  obtain  estimates  of  the  propagation  effects  on  optical 
beams  due  to  the  turbulence  in  a  fire  plume.  We  consider  the  propaga¬ 
tion  geometry  shown  in  Figure  36.  A  laser  beam  of  initial  diameter  D 
passes  through  a  plume  whose  centerline  is  halfway  between  the  laser 
and  the  screen.  As  before  the  propagation  path  is  horizontal  and  the 
point  of  closest  approach  of  the  beam  centerline  to  the  plume  center- 
line  is  specified  by  the  impact  parameter  distance  d.  On  that  part  of 
the  propagation  path  outside  the  plume,  the  beam  is  propagating  in  the 
ambient  atmospheric  turbulence.  As  shown  in  the  previous  subsection 
the  strength  of  the  index  of  refraction  fluctuations  in  the  ambient 
atmosphere  fs  expected  to  be  several  orders  of  magnitude  below  the 
plume  fluctuations.  For  relatively  short  path  lengths  1000  m  or 
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less)  we  expect  the  plume  turbulence  effects  to  completely  dominate  the 
atmospheric  turbulence  effects. 

The  plume  turbulent  parameters  required  for  propagation  calculations 
are  the  magnitude  and  spectrum  of  the  index  of  refraction  fluctuations 
is  a  measure  of  the  magnitude.  The  spectrum  depends  upon  the  char 
acteristics  of  the  plume  flow.  In  fully  developed  turbulence  of  the 
dissipative  type,  the  spectrum  of  the  velocity  fluctuations  is  known  to 
be  Kolmogorov.  The  ambient  atmospheric  turbulence  is  generally 
Kolmogorov.  A  conservative  passive  additive  is  a  quantity  which  is 
carried  along  by  the  turbulent  velocity  flow  tields  and  which  qocs  not 
influence  the  flow.  The  additive  fluctuations  assume  the  same  -  ■  ■ 
as  the  velocity  spectrum.  In  the  ambient  atmosphere  the  pour,:  ...1 
temperature  is  a  conservative  passive  additive.  Thus  the  ambient 
atmospheric  temperature  fluctuations  are  also  Kolmogorov.  Sin_e  in 
the  ambient  atmosphere  the  index  of  refraction  fluctuations  are  pro¬ 
portional  to  the  temperature  fluctuations  (Equation  51),  the  norma] 
atmospheric  index  of  refraction  fluctuations  are  also  Kolmogorov. 

In  a  fire  plume  the  initial  flow  is  convective,  decaying  into  a 
dissipative  flow.  The  velocity  fluctuations  in  the  dissipative  regime 
are  expected  to  be  essentially  Kolmogorov.  Although  the  plume  tempera¬ 
ture  is  not  a  passive  additive  since  the  plume  flow  is  driven  by  the 
buoyancy,  the  temperature  fluctuations  are  also  expected  to  be  essentially 
Kolmogorov  in  the  dissipative  regime.  As  we  have  seen,  the  index  of 
refraction  fluctuations  at  high  temperatures  are  not  directly  propor¬ 
tional  to  the  temperature  fluctuations  (see  Equation  53).  Thus  the 
index  of  refraction  fluctuation  spectrum  will  be  a  somewhat  modified 
version  of  the  temperature  spectrum.  However  since  we  do  not  know  the 
temperature  spectrum,  we  will  make  the  simplifying  assumption  that  in 
the  dissipative  flow  regime  of  the  plume  the  index  of  refraction  fluct¬ 
uations  have  a  Kolmogorov  spectrum.  Since  we  also  do  not  know  the 
spectrum  in  the  convective  region,  we  will  make  the  further  simplifying 
assumption  that  the  spectrum  is  Kolmogorov  everywhere.  This  assumption 
simplifies  the  propagation  calculations.  We  will  make  estimates  of  the 
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magnitude  of  the  fluctuations  depending  on  whether  the  flow  is  con¬ 
vective  or  dissipative,  but  we  will  use  the  same  spectral  shape 
(Kolmogorov)  in  both  flow  regimes. 

There  are  three  characteristic  size  regimes  associated  with  the 
turbulent  eddies.  According  to  the  Kolmogorov  theory,  these  regimes 
can  be  characterized  by  two  size  scales:  the  outer  scale  of  turbulence, 
Lq,  and  the  inner  scale  (also  called  microscale)  of  turbulence,  lQ. 

The  first  size  regime  is  called  the  input  range  and  here 

eddy  size  >  Lq 

It  is  in  this  regime  that  the  energy  is  extracted  from  the  mean  flows 
and  converted  to  turbulent  energy.  The  spectrum  in  this  range  de¬ 
pends  on  the  particular  flow  characteristics  and  no  general  formulas 
exist.  The  spectrum  is  normally  anisotropic  in  the  input  range.  The 
second  size  regime  is  called  the  inertial  subrange,  and  here 

(L  <  eddy  size  <  L 
o  o 

This  is  the  eddy  cascade  range.  The  large  eddies  break  up  into  smaller 

eddies,  which  in  turn  further  subdivide  into  still  smaller  eddies,  and 

so  on.  As  the  subdivision  continues,  any  initial  anisotropy  is  lost 

and  the  turbulence  is  essentially  isotropic  in  this  range.  The  inertial 

subrange  is  essentially  energy  conserving;  the  energy  loss  due  to 

viscosity  dissipation  is  negligible  compared  to  the  kinetic  energy  of 

an  eddy.  Kolmogorov  considered  the  physical  characteristics  of  fully 

developed  turbulence  in  the  inertial  subrange  and  deduced  that  the 

-11/3 

spectrum  was  proportional  to  K  where 

K  =  eddy  wavenumber  (m "*) 

-  2tt 
eddy  size 

The  third  size  regime  is  the  dissipation  range  and  here 
eddy  size  <  l 
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When  the  subdividing  eddies  reach  a  size  of  the  smallness  of  the  inner 
scale,  the  dissipation  of  the  eddy  kinetic  energy  due  to  viscosity  ef¬ 
fects  begins  to  dominate.  The  eddy  no  longer  subdivides  and  ceases  to 
exist  as  an  entity  as  the  viscosity  converts  the  eddy  kinetic  energy 
into  random  heat  energy.  The  spectrum  rapidly  drops  to  zero  in  this 
size  range. 

For  mathematical  convenience  the  three  regimes  are  often  combined 
in  one  spectral  formula,  usually  called  the  modified  von  Karman  spectrum 

VK>  ■  °-033  cn(K2*  tt)  "/6  p 

o'  '  m 

where 

<t>  =  spectrum  of  index  of  refraction  fluctuations 

2 

C  =  structure  constant  of  index  of  refraction  fluctuations 

"  c»-3'3, 

K.  ’  S.92/to  On’1). 

Figure  37  shows  a  sketch  of  the  spectrum.  The  input  range  is  shown 
dotted;  althouph  Equation  71  mathematically  covers  the  input  range, 
the  expression  must  be  considered  only  an  approximation  since  the 
actual  spectrum  depends  on  the  flow  characteristics  and  is  generally 
unknown . 


The  strength  of  the  turbulence  is  specified  by  the  magnitude  of 

2  2 

the  structure  constant,  C  .  In  the  ambient  atmosphere  C  varies  from 

-18  -2/3  ^  -13  -2/3  ^ 

about  10  m  for  weak  turbulence  to  10  m  '  for  strong 

turbulence.  For  the  Kolmogorov  spectrum,  the  relation  between  the 

structure  constant  and  the  previously  considered  index  of  refraction 

variance,  a2,  is  (References  13  and  14) 
n  . 

,  1.91  a 

r2  _  ft  2/3 

Cn  .2/3 “  m  •  (/2) 


o 

We  will  use  two  different  methods  to  compute  the  structure  constant 
in  the  fire  plume.  In  the  first  method  we  utilize  the  experimental 
measurements  of  Reference  10  and  assume 


80 


1/  _  _ 2tt _  i 

eddy  size  0 

Figure  37.  Kolmogorov  index  of  refraction  fluctuation  spectrum. 


oT(r)  =  0.4(T(r)  -Ta)  K  .  (73) 

2 

We  then  use  the  results  }f  the  previous  subsection  to  calculate  a  (r) 

2  n 

and  then  use  Equation  72  to  calculate  C^fr).  This  method  is  appropriate 
to  well  .mixed,  well  developed  turbulence.  We  would  expect  this  formal¬ 
ism  to  apply  after  the  initial  convective  plume  flow  has  decayed  away 
and  the  flow  has  become  dissipative  (ie,  Kolmogorov).  This  should  ap¬ 
ply  in  the  higher  altitude,  weakly  buoyant  regions  of  the  plume,  pro¬ 
vided  there  are  no  intermittency  effects  due  to  ambient  turbulence 
interactions. 
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In  the  second  method  of  calculating  C“,  we  relate  the  strength  of 
the  turbulence  to  the  gradients  of  the  mean  plume  parameters.  Assumiu 
the  temperature  is  also  essentially  Kolmogorov  in  spectrum,  we  have  an 
equation  analogous  to  72  for  temperature 

2  K91aT  2  -2/3 

C7,  =  - K  m  '  •■"4 

\J~  ^ 
o 

2 

where  is  the  structure  constant  for  the  temperature.  The  f  emperat 
structure  constant  for  Kolmogorov  fluctuations  can  be  re  l  at  r  !  •  rh> 
mean  temperature  gradients  by  (Reference^  13  an!  11 

C“  =  b  L  ' J  | grad  T|“  K  m  J  '7 

where  T  is  the  mean  temperature  and  b  is  a  constant  wh:Ji  <>  • 

the  flow  type.  The  exact  magnitude  of  b  is  not  we  <  !  1  now, .  .•> 
ambient  atmospheric  flow  h  is  estimated  to  be  in  the  ranee  ' 

3.5  (References  14  and  15).  Reference  IS  recommends  a  value  of  2."-. 
For  non-well  mixed  plume  flows  we  adopt  the  value 


b  =  3 


In  this  gradient  method  we  will  also  calculate  the  structure  constant 
using  a  second,  limiting  value  of  b  of 

b2  =  0.3  .  ("7) 

We  expect  this  second  value  to  be  appropriate  to  well  mixed  flows.  If 


we  approximate 


,  T  -  T 

2 


|  grad  T |  » 


and  assume  for  well  m  .  flows 


aT  *  0.4(T-Ta) 

wh<-'e  T  is  an  average  mean  temperature  over  the  plume  cross  section, 
''hen  using  Equations  74  and  75  we  find 
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b2  «  1.91  x  (0.4)2  =  0.3 

Combining  Equations  74  and  75  we  have  in  the  general  case 
?  b  L2|grad  T(r,z)|2 

a2(r,z)  =  - 2 - - -  K  .  (78) 

For  the  plume 

igradT(r>I)i^(^i)2.(m^)2  (7., 

and  the  radial  and  vertical  gradients  at  a  point  are  found  from  the 
plume  model  of  Section  2.  Lq  is  the  outer  scale  of  plume  turbulence, 
the  length  of  the  temperature  gradient  region.  For  plumes 


L  «  R  m 
o 


(80) 


In  this  second  method  we  assume  non-well  mixed  flows  and  use  Equations 

2  2 

69  and  70  to  convert  from  the  of  Equation  78  to  a^.  Equation  72 
then  gives  C2. 

2 

With  Cn  defined,  we  next  develop  the  propagation  formulas.  A 
primary  measure  of  turbulence  effects  is  the  coherence  length.  The 
coherence  length  of  a  spherical  wave  after  propagation  from  the  screen 
to  the  laser  is  (References  15  and  16) 


=  [‘,e  A.j 


d£(l  -  C)5/3  C2(£xs 


■3/5 


(81) 


where 


Pq  =  coherence  length  (m) 
x  =  distance  from  laser  (m) 
xg  =  total  distance  from  laser  to  screen  (m) 
k  =  2ir/Am  (nf1) 

X  =  wavelength  of  laser  (m) 

C  (x)  =  structure  constant  of  index  of  refraction  fluctuations  at 

n  -2/3 

point  x  (m  ) . 
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We  assume  that  in  the  ambient  atmosphere  the  turbulence  is  constant 
and  equal  to 


2  -14  -2/3 

(C  1 ,  =  constant  =10  m 
n  A 


a  value  representative  of  relatively  strong  ambient  turbulence.  The 
ambient  turbulence  effects  are  generally  negligible  compared  to  the 
plume  turbulence  effects  for  relatively  short  path  lengths  (<1000  r; s 
We  have  written  a  computer  program  to  solve  Equation  SI  ior  the 
Gaussian  plume  of  Section  2.  Far  from  the  plume  C~(x)  ihe  atmo¬ 
spheric  value;  inside  C“  (x)  is  the  plume  value;  at  the  fringes 
2  D 

plume  we  take  C^(x)  as  the  maximum  of  the  atmospheric  or  plume  value*-- 

Next  consider  the  size  and  location  of  the  beam  spot  site  ;n  the 
screen.  The  effect  of  turbulence  will  cause  the  beam  to  spread  ;in 
addition  to  the  vacuum  spread]  and  to  move  about  from  place  to  pi  act¬ 
on  the  screen  (to  "jitter").  Figure  38  shows  a  schematic  of  the  spot 
geometry  at  the  screen.  At  the  screen  the  instantaneous  spot,  of 
average  radius  ps,  jitters  from  place  to  place  on  the  screen  with  a 
time  constant  of  the  order  of  D/v,  where  D  is  the  beam  diameter  and 
v  is  the  rise  velocity  of  the  plume.  The  average  distance  of  the 
centroid  of  the  instantaneous  spot  from  the  center  of  the  screen  is 
Pc>  For  a  long  exposure  time  (average  over  many  locations  of  the 
instantaneous  spot)  the  average  radius  of  the  smeared  spot  is  p^. 


First  consider  the  propagation  of  the  laser  beam  in  free  space. 

At  the  aperture  (x = 0)  we  assume  the  amplitude  distribution  is  Gaussian 
with  a  beam  size  (radius)  of  D/2  and  the  phase  distribution  is  parabolic 


with  a  radius  of  curvature  -F  .  This  phase  distribution  will  produce  a 


beam  focused  at  x=F 


distance  x  from  the  laser  is 

s  .  2  i 

,  4x  „2  /  x 


The  free  space  beam  radius,  p^,  at  the  screen 

2 


CD  = 


s  D"  /.  s\ 

■D2  +  4  '  V 


(82, 


This  is  the  diffraction  limited  behavior  of  a  perfect  laser  beam.  In 


our  calculations  we  consider  a  collimated  beam  for  which  F  =  °°. 
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Figure  38.  Spot  geometry  at  screen. 

The  mean  square  of  the  long-term  beam  radius  is 
(Reference  16  ) 


CD 

// 


d2p  p2r2(*s,p,p) 


// 


d  p  r2(xs,p,p) 


2 

m 


p  =  radial  vector  lying  in  the  plane 


usually  defined  by 


(83) 


of  the  screen  (m) 


rfhere 


—+■  •+ 

T.,  (  x,  Dj  ,p.,)  =  second  moment  of  the  electric  field 
=  <u(x,Pj)u*(x,p2))y 

u  ~  electric  field  amplitude 
^indicates  ensemble  average. 

is  commonly  called  the  mutual  coherence  function.  For  the  Ko; 
spectrum  (Reference  lb) 

.  4x2 


'  PD  +  ,  2  2 


<PL>  = 


k  P, 


6.6xJ 

*  p/3  / 

o  J  0 


dC Cl  -  02Cn2(C*s 


where 


\  =  [k2C2l5/V!  m 

L  n  o 

£  =  inner  scale  of  turbulent  fluctuations  tm) . 

o 

In  the  atmosphere  iQ  is  of  the  order  of  a  millimeter  (10  m.  and  we 
assume  this  value  for  our  calculations.  The  actual  value  o* 
pends  upon  the  rate  at  which  the  turbulent  energy  is  being  generated 
and  dissipated.  The  computer  routine  which  evaluates  Lquation 
81,  by  numerical  integration  also  evaluates  the  integral  in  the  sec 
ond  part  of  Equation  84. 


The  mean  square  centroid  location  can  be  written 


d  a,  J  4(xs.0rpj  ,P2,P2) 


(85) 


// 


dVi  r:(W:V 


where  T  is  the  fourth  moment  of  the  field 
4 

r4(x  ,0  ,p2,03,p4)  =  <^u(x  ,P  )u*  (x  ,P2)u(x  ,P3)U  (  x , P4 

With  ar>d  /p^^>  defined,  the  definition  of  <(p2^>  fol  1  ows  from  the 

relation 
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(86) 


<?2i>  ■  <p2>  *  <?5> 

A  rigorous  solution  exists  for  the  second  moment  V ^ ,  so  <(p^  can 
be  calculated  directly.  Unfortunately  there  is  no  rigorous  solution 
for  the  fourth  moment  so  ^p2^>  and  cannot  be  calculated  direct¬ 

ly.  However,  there  are  a  number  of  approximations  available  for  F,, 
/p^>  ,  and  <PS^>  •  We  adopt  the  approximations  given  by  Fante  in  Ref¬ 
erence  16.  Fante  presents  a  number  of  separate  cases  depending  on  the 
strength  of  the  turbulence  and  the  path  length.  Define 

2  2 

x^  =  minimum[kp^,kD  ]  m 

Case  1 : 


p  <<  D  <  L  and  x  <  x  _ 
o  o  s  ~  pD 


Then  approximately 
S.2' 


2  4xs 

/P 

l/3-i 

s  PD  +  .2  2 
k  Po 

1  -  0.62U 

r)  J 

2.97xJ 

2 

m 

6/5 


k  pJ/ JD 
o 


For  Case  1  we  use  Equation  86  to  calculate  ^>2^ 
Case  2: 


m2  (87) 

(88) 

instead  of  Equation  84. 


p  ~  D  and  x  <  x  n 
o  s  ~  pD 

For  this  case  there  are  no  simple  analytic  expressions.  Define 


Then  <(ps/>  can  be  obtained  from  the  numerical  results  plotted  in  Fig¬ 
ure  39.  Knowing  and  <(p2^>  follows  immediately  from  Equation 

86. 
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Figure  39.  Ratio  of  the  short-  to  long-term-averaged 
beam  spread. 

Case  3: 


Po  »  D  and  xs  <  xpD 
For  this  case  there  is  very 

<Ps)  *  <0 


little  beam  jitter  and 


/  2\ 
<Pc> 

Case  4: 


0 


x  >>  x 
s  pD 


Again  for  this  case 


<pf>  -  <p2l> 


and  can  be  neKlected.  However  case  3  was  for  weak  turbulence. 

This  case  is  due  to  strong  turbulence.  Here  the  beam  has  lost  its  co¬ 
herence.  The  instantaneous  picture  on  the  screen  no  longer  consists  of 
a  single  spot  but  consists  of  a  multiple  number  of  smaller  string-like 
spots  at  random  locations  within  the  PL  radius. 


A  computer  routine  has  been  written  to  solve  these  propagation  equa¬ 
tions.  Results  of  some  sample  calculations  are  shown  in  Tables  3  through 
11  and  Figures  40  through  43.  The  tables  show  propagation  results  for  a 
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Table  3  La  ser  ; -trail  propagation  results  for  a  plume  with  initial  parameters 
T  -1  =  b00  k,  V  =10  m  s'1,  R=1  m. 
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Laser  beam  propagation  results  for  a  plume  with  initial  parameters 
T„-T  =1000  K,  V  =20  m  s"1,  R=10  m. 
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Laser  bean  propagation  results  for  a  plume  with  initial  parameters 
T  -T  =1500  K,  V  =20  m  s'1,  R=10  m. 
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COHERENCE  LENGTH  (m) 


two  meters  above  the  plume  reference  plane  (flame  tips)  which  intersects 
the  plume  centerline.  Initial  plume  radii  were  1,  5,  and  10  meters  with 
initial  rise  velocities  of  10,  15,  and  20  m  s  *,  respectively.  For  each 
plume,  three  initial  excess  centerline  temperatures  of  500,  1000,  and 
1500  K  were  taken.  These  cases  correspond  to  the  plumes  shown  in  Figure 
4,  7,  10,  13,  16,  19,  22,  25,  and  28.  Laser  wavelengths  were  chosen  at 
1  and  10.6  pm;  the  beam  was  taken  as  collimated  with  initial  beam  diam¬ 
eters  of  1  mm,  1  cm,  and  10  cm.  The  two-order  magnitude  variation  in 
beam  diameter  was  chosen  for  two  reasons: 

1.  to  illustrate  the  sensitivity  of  the  results  to  beam 
diameter,  and 

2.  to  cover  the  expected  diameter  range  of  laboratory, 
experimental,  and  battlefield  lasers. 

Total  path  lengths  chosen  were  30,  100,  and  1000  meters,  with  the  plume 
centered  at  the  midpoint  of  the  path.  The  path  length  range  was  chosen 
to  cover  the  typical  shorter  paths  of  an  experimental  test  and  the  longe 
paths  of  the  battlefield. 

Results  are  given  for  the  three  mixing  models  where  for 
Mixing  model  1 

2 

3:  b  =  1 

Mixing  model  1  is  the  well -mixed,  Gaussian  fluctuation  case.  Mixing 
models  2  and  3  are  the  non-well  mixed  .l^crete  bimodal  temperature 
fluctuations.  The  turbulent  propagation  results  tabulated  are  the 
coherence  length,  and  root  mean  square  values  for  long-  and  short-term 
beam  spread,  and  beam  centroid  offset.  Also  given  for  con^iarison  pur¬ 
poses  is  the  vacuum  spread  of  the  beam,  p^;  the  difference  between  the 
total  and  vacuum  spread  is  due  to  turbulence  effects. 


'J_  =  0.4  (T  T  ) 

*  u 

b  I.*" 1  grad  I 
_»  o 

T  1.91* 


b  =  0.3 
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Figures  40  through  43  illustrate  the  variation  with  altitude  and 
radial  offset.  For  the  plume  of  Figure  16,  initial  conditions  Tc~Ta  = 
1000  K,  V  =  15  m  s'1,  R  =  5  m,  the  propagation  parameters  were  cal¬ 
culated  as  a  function  of  altitude  for  the  100-m  path  length  for  the 

1- prn  wavelength,  and  the  1-cm  beam  diameter.  All  paths  intersect  the 
plume  centerline.  Results  are  shown  in  Figures  40  and  41  for  mixing 
models  1  and  3.  To  illustrate  the  radial  offset  effect,  the  path  at 

2- meters  altitude  was  progressively  offset  from  the  plume  centerline. 

The  offset  results  are  shown  in  Figures  42  and  43,  again  for  mixing 
models  1  and  3.  The  dotted  segments  in  Figure  40  is  the  region  of 
transition  from  one  propagation  approximation  to  the  next  (from  case  1 
to  case  2,  Equations  87  to  90).  The  long-  and  short-term  radii  are 
relatively  continuous,  but  the  centroid  offset  is  slightly  discontinuous. 
Over  the  first  25  meters  of  this  plume  there  is  a  relatively  small 
change  in  the  propagation  parameters  with  altitude.  The  excess  center- 
line  temperature  has  fallen  from  1000  K  at  the  reference  plane  to  only 
about  200  K  at  25  meters.  As  shown  in  Figure  35,  200  K  is  still  near 
the  peak  of  the  index  of  refraction  fluctuations.  The  turbulence  ef¬ 
fects  will  diminish  rapidly  with  temperature  as  the  temperature  drops 
below  about  100  K.  The  turbulent  effects  for  the  non-well  mixed  flows, 
mixing  model  3,  are  about  a  factor  of  3  larger  than  the  well-mixed 
Gaussian  model  1.  This  is  about  the  range  of  effects  to  be  expected 
between  the  convective  and  dissipative  regions  of  the  plume. 

As  expected  the  turbulence  effects  decrease  rapidly  as  the  ray 
path  approaches  the  plume  edge  (Figures  42  and  43).  Not  only  is  the 
plume  temperature  decreasing,  but  the  total  path  length  through  the 
plume  turbulence  is  also  decreasing.  Again  the  difference  in  mixing 
models  is  about  a  factor  of  3  in  propagation  effects. 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  this  report  we  have  developed  models  for 

1.  the  mean  plume  parameters  including  radius,  temperature, 
and  velocity  as  a  function  of  altitude, 

2.  the  turbulence  associated  with  the  plume  mean  parameters, 
and 

3.  the  effects  on  optical  beams  propagating  through  the  plumes 
including  both  mean  bending  and  turbulence  effects. 

Some  of  the  models  have  been  verified  by  previous  experimental  evidence, 
but  verification  of  the  remaining  models  must  await  the  planned  plume 
test  series. 

The  models  for  the  mean  radius,  temperature,  and  velocity  of  the 
plume  as  a  function  of  altitude  have  a  solid  theoretical  base,  and 
have  been  verified  by  experimental  measurements.  We  expect  the  results 
predicted  by  the  mean  plume  models  to  be  fairly  accurate  when  compared 
to  experimental  data  on  essentially  ideal  plumes.  An  ideal  plume  is  a 
plume  generated  under  controlled  conditions  with  negligible  intermittency 
effects.  An  ideal  plume  can  be  readily  generated  under  laboratory  con¬ 
ditions,  but  may  be  difficult  to  generate  under  field  conditions.  As 
mentioned  previously,  interactions  of  the  plume  flow  fields  with  the 
ambient  flow  fields  may  introduce  significant  intermittency.  The  plume 
may  move  about  in  space  or  may  even  break  up.  If  the  plume  merely  moves 
about,  then  the  model  predictions  can  be  used,  noting  that  the  predic¬ 
tions  apply  at  the  instantaneous  location  of  the  plume.  But  if  the 
plume  breaks  up  or  the  ambient  interactions  introduce  significant  addi¬ 
tional  mixing,  then  the  model  predictions  would  no  longer  apply. 
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In  contrast  to  the  fairly  solidly  based  models  for  the  mean  plume 
parameters,  the  models  for  the  plume  temperature  fluctuations  are 
primarily  theoretical  with  little  experimental  verification.  For 
well  developed  plume  turbulence,  experimental  data  verify  the  well 
mixed  formulation  model,  Equation  61;  see,  for  example.  Reference  10. 

But  the  non-well  mixed  models  are  purely  theoretical  with  no  experimental 
verification.  Moreover,  for  both  the  well  and  the  non-well  mixed 
models,  there  is  no  experimental  verification  for  the  assumption  that 
the  index  of  refraction  fluctuations  have  a  Kolmogorov  spectrum.  We 
can  reasonably  expect  the  well  developed,  well  mixed  flow  to  be 
Kolmogorov.  But  until  the  initial  convective  flow  decays  into  a  dis¬ 
sipative  flow,  we  do  not  expect  the  spectrum  to  be  Kolmogorov.  Hope¬ 
fully  the  planned  test  series  will  measure  the  spectrum. 

For  small  total  bending  angles,  the  ray  tracing  models  for  the  mean 
bending  are  solidly  based  in  both  theory  and  experiment.  The  problem 
lies  in  defining  the  mean  thermal  environment  of  the  plume;  the  tracing 
of  a  ray  through  this  environment  is  straightforward. 

Given  that  the  turbulence  spectrum  is  Kolmogorov,  the  propagation 
models  have  been  derived  theoretically  and  verified  experimentally  (see 
Reference  16).  However  the  experimental  verification  has  been  primarily 
for  propagation  in  ambient  atmospheric  turbulence.  We  also  expect  the 
models  to  apply  to  propagation  through  the  localized  highly  intense 
plume  turbulence,  but  experimental  verification  awaits  the  upcoming 
test  series. 

An  experimental  test  series  on  fire  plumes  will  be  conducted  at  the 
White  Sands  Missile  Range  in  the  near  future.  This  test  series  should 
eliminate  many  of  the  current  uncertainties  in  the  plume  models.  The 
present  plume  and  propagation  models  can  be  used  for  test  planning. 

For  field  test  propagation  predictions,  we  recommend  mixing  model  3, 
the  non-well  mixed  flow  model.  The  initial  plume  flow  is  convective, 
which  is  non-well  mixed.  As  the  altitude  increases  the  convective  flow 
decays  into  a  dissipative  flow.  If  fully  established  with  no  inter- 
mittancy,  the  dissipative  flow  should  be  well  mixed,  corresponding  to 
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mixing  model  1.  But  for  the  field  tests,  we  expect  the  ambient  inter¬ 
actions  to  introduce  intermittances,  tending  to  produce  non-well  mixed 
flows.  Thus  for  field  tests,  mixing  model  3  is  recommended  for  all 
altitudes  for  planning  purposes. 

The  specific,  details  of  the  forthcoming  plume  test  series  will  be 
determined  by  the  available  test  equipment  and  budget  considerations, 
but  some  general  recommendations  can  be  made.  A  variety  of  fire  sites 
and  burning  rates  should  be  employed  if  time  and  budget  allow.  However 
it  is  more  important  to  take  data  at  a  number  of  altitudes  in  a  single 
plume  than  to  take  data  at  a  single  altitude  in  a  number  of  different 
plumes.  To  eliminate  the  complicating  effects  of  particulates,  the 
fires  should  be  clean,  i.e.,  essentially  smokeless. 

Ideally,  over  the  data  collection  period  for  a  given  plume,  the 
intrinsic  parameters  of  the  plume  should  be  time  independent .  The 
burning  rate  of  the  fire  should  be  constant.  The  ambient  conditions 
should  be  chosen  as  calm  as  possible  to  minimize  ambient  interactions. 

If  there  is  significant  lateral  motion  of  the  plume,  then  measurement 
should  be  recorded  relative  to  the  instantaneous  plume  centerline. 

This  relative  measurement  requires  that  the  probe  or  laser  beam  posi¬ 
tion  and  the  plume  position  be  recorded  at  each  measurement  time. 

To  characterize  the  mean  plume  parameters,  the  plume  diameter  and 
mean  centerline  temperature  and  rise  velocity  should  be  measured  at  a 
number  of  altitudes.  At  selected  altitudes  the  temperature  and  velocity 
should  be  measured  as  a  function  of  horizontal  distance  from  the  plume 
centerline. 

Since  the  ultimate  purpose  of  the  plume  modeling  effort  is  to  pre¬ 
dict  (. "opagation  effects,  measurements  of  various  propagation  effects 
are  a  must.  For  direct  comparison  with  model  predictions,  measurements 
should  include  short- and  long-term  beam  spread  and  jitter  (centroid 
location).  Measurement  of  the  short-term  beam  spread  and  jitter  re¬ 
quires  that  the  response  time  of  the  measuring  system  be  fast  enough  to 
measure  the  instantaneous  beam  spot  as  it  flits  from  point  to  point  on 
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the  receiving  screen.  The  time  it  takes  the  spot  to  move  from  one  point 
to  the  next  is  approximately 


v 


where  D  is  the  diameter  of  the  beam  at  the  plume  location  and  v  is  the 
plume  rise  velocity. 

The  diameter  of  the  beam  when  it  passes  through  the  plume  should  be 
small  compared  to  the  diameter  of  the  plume.  This  will  minimize  the 
effects  of  inhomogenities  in  the  turbulence  strength  and  spectrum  over 
the  cross  sectional  dimension  of  the  beam.  Optical  beams  should  have 
no  problem  meeting  this  smallness  criterion.  However  if  microwave 
beams  are  fielded  as  an  adjunct  to  the  optical  tests,  some  care  may  be 
required  in  the  relative  placement  of  transmitter,  plume,  and  screen 
(receiver) . 

The  temperature  is  the  controlling  system  parameter  for  plume 
turbulence  and  propagation  effects.  At  optical  frequencies  the  tempera¬ 
ture  fluctuations  determine  the  index  of  refraction  fluctuations;  the 
contributions  due  to  pressure  and  water  vapor  fluctuations  are  negligible. 
The  quantity  required  for  turbulence  propagation  calculations  is  the 
spectrum  of  the  index  of  refraction  fluctuations,  which  in  turn  follows 
directly  from  the  spectrum  of  temperature  fluctuations.  If  the  instru¬ 
mentation  is  available,  measurements  of  the  plume  temperature  spectrum 
would  be  extremely  valuable.  A  series  of  measurements  at  increasing 
altitudes  would  determine  the  decay  of  the  initial  convective  flow  into 
the  dissipative  flow.  The  measurements  will  determine  if  or  when  the 
index  of  refraction  fluctuations  approach  the  presently  assumed  Kolmogorov 
spectrum. 
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